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THE MIXING OF A JET OF GAS WITH AN 
ATMOSPHERE OF A DIFFERENT GAS AT 
LARGE DISTANCES FROM THE ORIFICE 


PART I. THE PLANE JET 


By L. J. CRANE and D. C. PACK 


(Department of Mathematics, the Royal College of Science and 
Technology, Glasgow) 


[Received 31 October 1960} 


SUMMARY 

This paper is concerned with the flow of a plane jet of gas out of a narrow slit 
into a medium consisting of another gas. Approximate p. files of velocity and 
concentration in the jet are determined for both laminar and turbulent flows. The 
solution is valid only at large distances from the orifice. 

In order to solve this problem analytically four basic assumptions have been made, 
namely, that the gases are incompressible, i.e. that the velocity of the fluid is small 
compared with the local speed of sound, that the temperature is constant every- 
where, that the boundary layer approximations are valid, and that the pressure is 
constant over the mixing region. In the case of turbulent flow the Reichardt 
constant exchange coefficient hypothesis is used. 

In this paper a method similar to that used by Pack (3) for the axially symmetric 
compressible jet is employed. This method as applied here consists in expanding 
both the stream function and the concentration in a Rayleigh—Janzen series of 
powers of a parameter whose magnitude depends on the concentration of one gas 
in the mixture. 


1. Introduction 


THE first study of the mixing of jets of inhomogeneous fluid was made by 
Chou (1) in 1947. Chou considered the equation of diffusion along with 
the boundary layer equations. Chou’s work was extended by Pai (2), who 
transformed the diffusion and boundary layer equations into equations 
of heat-conduction type. Pai solved these equations numerically for flow 
in which the jet was embedded in a stream moving with a uniform velocity 
not very different from that of the jet. 

This paper is concerned with the flow of a plane jet of gas out of a narrow 
slit into a medium consisting of another gas. Approximate profiles of 
velocity and concentration in the jet are determined for both laminar and 
turbulent flows. The solution is valid only at large distances from the 
orifice. 


(Quart. Journ. Mech. and Applied Math., Vol. XIV, Pt. 4, 1961] 
5092 .56 ce 
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In order to solve this problem analytically four basic assumptions have 
been made, namely, (i) that the gases are incompressible, i.e. that the 
velocity of the fluid is small compared with the local speed of sound; 
(ii) that the temperature is constant everywhere; (iii) that the boundary 
layer approximations are valid; and (iv) that the pressure is constant over 
the mixing region. In the case of turbulent flow the Reichardt constant 
exchange coefficient hypothesis is used. 

In this paper a method similar to that used by Pack (3) for the axially 
symmetric compressible jet is employed. This method as applied here 
consists in expanding both the stream function and the concentration in 
a Rayleigh—Janzen series of powers of a parameter whose magnitude 
depends on the concentration of one gas in the mixture. 


2. Equations of motion 


Let u and v be the velocity components (or mean velocity components 
in the case of turbulent flow) parallel to rectangular cartesian axes O(z, y). 
Let the origin of coordinates be taken on the axis of the jet near to the 
point at which mixing begins, i.e. near the orifice. It can be shown that, 
if the pressure and temperature are constant in the mixing region, then 
the density p of the gas mixture in the jet is given by 


p* = p/p, = 1+8e, 
where p, is the density of the gas in the surrounding medium (p, is the 
density of the gas 1 coming out of the orifice); 8B = (m,/m,)—1, m, and 
m, being the molecular weights of the two gases, respectively; and c is 
the molecular concentration of gas 1 in the mixture. When c is small the 
viscosity u of the mixture of gases in the jet is given by 

p* = p/p, = 1+y¢, 
where y = (,/A,p2)—A,, , and p, are the viscosities of gases 1, 2 
respectively and A,, A, are constants of the gases [cf. (4)]. 

The equations of motion for both laminar and turbulent flow of a steady 

plane jet whose main direction of motion is that of x increasing, are (5) 


not pou é ou (1) 
— j—- == —- |] €9 — 
at ae Cy oy Poy 


e é 
and — (9us)-+- (pe) == 0, 9 
5 Pu) dy (P ) ) 


In laminar flow « = /p, while in turbulent flow « is an empirically deter- 
mined average value of the eddy kinematic viscosity over a section of the 
jet, which depends only on z. 

In addition to (1), (2) it is also necessary to take into account the 
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equation governing the diffusion of the gas in the jet into the surroundings. 
This equation is 


é é al ec “ 
=~ (cu)+ — (ev) = — ‘ox; (3) 
Cx oy oy oy 


where ¢, is the molecular coefficient of diffusion, D,. in the case of 
laminar motion; D,, is taken here to be constant following Jeans (4); while 
in the case of turbulent flow ¢, is an empirically determined mass transfer 
coefficient defined in a similar way to €, €, being assumed to be a constant 
multiple of «; the constant of proportionality Z is defined by «ep = Ee. 

The equation of continuity, (2), can be satisfied by introducing a stream 
function % defined by 

pu=y,, pu= —yp,. 
Equations (1) and (3) are simplified by changing the independent variables 
from (x, y) to (C,z) where - 
t= ( p* dy, 
0 


and ¢ stands for x when the flow is laminar, while if the flow is turbulent 


a 


. {)q 
dx 


€o 


where ¢, is Reichardt’s constant exchange coefficient. In terms of these 
variables (1) now becomes 


eh Oh ey Oh _ 2 [41024] 


dz Oldz aC a2? ez? 


(4) 


where d(c) = u*p* in laminar flow and ¢(c) = (p*)* in turbulent flow. 
é(c) is approximately equal to 1+-nc when c is small. In laminar flow 
¥ = p, and nm = B+; when the flow is turbulent « = €,p, and n = 28. 
The equation of diffusion now becomes 
fp Oc op Oc éc ‘ 
_ > = = ao(p*)? — S|e*s | (5) 
Oz a al éz ez 
where o is a dimensionless constant of order unity which is equal to 
Dy» py) 4g, the Schmidt number, when the flow is laminar and E when the 
flow is turbulent. 


3. Boundary conditions 
The boundary conditions on the velocity and concentration are: 


on the axis of the jet, y = 0: 


oa & mu  & 9 


ey ey 


> 
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and on the boundaries of the jet, y = +00: 
u= 0, } az @, 
In terms of the new variables these become: 


éc é 
onz=0: =0, —=0, and —=0, 
62? 


éz 
a 
andonz=-+0: c=0 and est 
éz 


4. Solutions of (4) and (5) 

A first approximation to the solution of (4) is obtained by putting p* = 1, 
i.e. by neglecting variations in density. When this approximation is made 
the stream function is (6) 

6aaf' tanh é, £ = az/{i, (6) 
where a is a constant of the motion which depends on the initial momen- 
tum flux. 

A fir: : approximation to c may now be obtained by inserting the 
approximate stream function (6) in (5) and neglecting the terms of second 
order and higher in ¢ in that equation. Equation (5) then becomes 

e a a a 22 
+ thon -_ 10. (7) 
dz Of = a a2? 
A solution of (7) of the form 
] 
om 9(€) 
a 
is now sought. The value of m is found from the condition that the flux 
of molecules of gas 1 crossing any section of the jet is a constant which 
is independent of x. This requirement determines that m = }, i.e. 


¢= 


c= an). 


When a concentration function of this type is inserted in (7) the required 


solution is b 
c= Zi (sech2é)"2, 


where 6 is a constant whose magnitude depends on the flux of molecules 
of gas 1 issuing from the orifice. Its value is calculated using the result 
that the total mass flux, M, of molecules of gas 1 from the slit is 


M = 2(8-+1)p, { ewdy = Baab(B-+1)B(~ +1; i), 
0 


where B(/;m) is the beta function. 
It now remains to be seen what effect the diffusion of the first gas into 
the surrounding atmosphere produces on the velocity profile. The method 
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used is to expand the stream function and the concentration in a Rayleigh- 
Janzen series in powers of a parameter which should be related to the 
concentration in some way. In fact this parameter is the first approxima- 
tion to the concentration on the axis of the jet, namely, b{-. Thus the 
appropriate expansions for % and c are 

% = 6aal*{ tanh §+6f-*F(€)+...] (8) 
and c = b{-*{sech?°é+bf-4G,(€)+...], (9) 
use being made of the first-order approximations. When (8), (9) are in- 
serted in (4) and the terms of order b[-* equated to zero, F,(€) is found to 
satisfy the equation 


FY +2 tanh €F{+6sech*éF, = 2n _ltanh &(sech*£)!+"?], (10) 


where dashes denote derivatives with respect to ¢. Let ¢ = tanhé and 
P = dF,/dé, then (10) becomes after some simplification 
dP 
dt® 
An integral of the left-hand side of this equation is ¢(1—?*); thus the sub- 
stitution P = t(1—#*)H gives 

(1 —#)?H” +2(1—#)(1—38)H’ = 2n— [1 —#)'+¥], 


d 
]— 2 +6P = 2In—\ltH1— 2\l+Lo ' 
(1—f*) 5 na le {2)I+1e] 


This equation has the first integral 
#2(1—#2)2H’ = 2nt?(1—t?)!+Vo+ no 


+20 


(l—#)+¥° ©, 


where C is a constant whose value in this case is zero because of the con- 
dition that dF,/dé must vanish on the boundaries of the jet. Finally, the 
solution of (10) which satisfies the condition that Fj(é) = 0 on the axis 
€ = 0, is 
dF, 
dé 


oO 


= (1—#)H = —“*_ (1—#)[2tBy(4; 1/o)—(1—?)"*), 
1+26 


Zz 
where B.(p:q) | (L—u)*-twP-? du. 
0 
When o is unity the solution is simply 
dF, a 
dE 4n(1 —t?)(5t?—1). 


TABLE oF (1 —f*)(5t?—1) 





o"2 o3 o"4 os 





0°424 | 0°549 
0°168 | 0 188 
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5. Change of scale 

It can be shown that 

ayf-* = £, = €—bBl-* tanh é. (11) 
Thus the velocity profile in terms of £, is obtained by distorting the 
velocity profile in terms of € by a translation approximately equa! to 
—bBl- tanh é. 

6. Conclusion 

The x-component of the velocity is approximately 

6a*x b 


tl |sechte + zi Fi6| . 
2 


A first glance at this result shows that the fractional change in the axial 
velocity is approximately ( bi ve 


C1 +20) 

In laminar flow y is, in general, only about a third or a quarter of 8, so 
that the order of magnitude of the fractional change in velocity on the 
axis is about —4b8 (assuming o to be about unity). This result seems to 
be in disagreement with the physics of the problem. For it means that 
if the molecules of gas, issuing from the orifice, are heavier than the 
molecules of the surrounding atmosphere, then the velocity of the jet on 
the axis is less than it would be if both gases had the same density. It 
might have been expected that a jet of dense gas would have a sharper 
profile than a less dense gas, because a massive molecule will suffer smaller 
deflexions on average than a light molecule in random collisions with the 
molecules of a given gas. This paradox is resolved by remembering that 
the total mass flux and the total momentum flux have been kept constant. 
Thus when the gas issuing from the slit is made up of massive molecules 
its velocity is from the start reduced so that the total momentum flux may 
be kept constant. The velocity profile is therefore not really significant 
in this problem. 

A more meaningful characteristic of the motion is the momentum flux, 
pu®, which is approximately 
36a%a? sech*é ! + 6 

Pe ti vi 
when o = |. It is found, when the independent variable ¢ is changed to 
&, by means of (11), that the momentum flux profile in terms of the latter 
variable is not sensibly affected by differences in density in the flow, 
provided these are not too large. Take for example the case when 
Bbt-* = }. In this instance it is found that, for a jet of hydrogen going 


2n : ) 
[Bsech*¢-+ —g (Stanh*¢- i : 
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into an atmosphere of oxygen and vice versa, the variation in the momen- 
tum profile over the core of the jet is only 1 or 2 per cent. This result 
is true for both laminar and turbulent flows. (It is to be noted that when 
the flow is turbulent the solution predicts quite appreciable changes in 
the momentum profile near the edges of the jet. These, however, would 
be masked in practice by the intermittency of the turbulence on the edges 
of the jet.) Thus provided the density differences are not too great, say 
10 to 20 per cent at most, the momentum profile is sensibly that of Bickley, 
namely, an 
pue = ata et wachit€. 
pot! 

It is of interest that the conclusion that the momentum profile is inde- 
pendent of density differences has also been shown to be true for axially 
symmetric turbulent jets by the experiments, covering a wide range of 
gases and density differences, of Vulis and Terekhina (7). 

Finally, it is to be remembered that these results only hold when € is 
large, i.e. at large distances from the orifice, because it is only in this 
region that the original solution of Bickley, on which this work depends, 
is valid. 


“ork described here was submitted by L. J. Crane as part require- 


¢ the degree of Ph.D. of Glasgow University. It was carried out 
» was the holder of a Sir James Caird Travelling Scholarship. 
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SUMMARY 


The flow considered in this Part is that of an axially symmetric jet of one gas 
which diffuses into an atmosphere made up of a different gas. The solution is 
developed in an asymptotic series whose coefficients depend on the axial concentra- 
tion, in the gas mixture, of the gas which comes out of the orifice. 


1. Introduction 


Tus Part is concerned with the analogue, in axially symmetric flow, of 
the problem considered in Part I [Crane and Pack (1)]. That is to say 
the flow is that of an axially symmetric jet of one gas which diffuses into 
an atmosphere made up of a different gas. As in Part I, certain basic 
assumptions have been made so as to obtain the solution in an analytical 
form. The solution is developed in an asymptotic series whose coefficients 
depend on the axial concentration, in the gas mixture, of the gas which 
comes out of the orifice. 


2. Equations of motion 


Since the flow is of boundary layer type at constant pressure, the 
equations of motion are 


a ay 1 y aN 
of oes 4 a. a <(eor<) (1) 


Cx or ror or 


and Ps (pur)-+ es (pur) = 0. (2) 
Ox or 


A system of cylindrical coordinates (x,r) is chosen so that the z-axis lies 
along the forward axis of the jet; u and v are the componentst of velocity 
along and perpendicular to this axis respectively. The quantity «, in (1), is 
u/p in laminar flow and is an empirically determined coefficient of eddy 
kinematic viscosity when the flow is turbulent; » and p are the coefficient 
of viscosity and the densityt of the gas mixture, respectively. 

+ The mean value is to be taken in the case of turbulent flow. 


(Quart. Journ. Mech. and Applied Math., Vol. XIV, Pt. 4, 1961) 
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In addition to equations (1) and (2) it is necessary to take into account 
the diffusion of the first gas 1 (which issues from the orifice) into the 
atmosphere (which consists of gas 2). The boundary layer equation of 
diffusion is in this case 

é é é oc ‘ 

5, ert a (cur) = <leore}, (3) 
where ¢ is the molecular concentrationt of the gas 1 in the mixture. The 
quantity «, is the molecular coefficient of diffusion D,, in laminar flow 
and is an empirically determined mass transfer coefficient when the flow 
is turbulent. 

Experiments on round turbulent jets (2) show that at large distances 
from the orifice the non-dimensional velocity and concentration profiles 
are functions of r/x. When functions of this type are substituted in the 
equations of motion it is necessary for consistency that e and ¢, have the 
constant values «, and He, respectively. 

The continuity equation (2) is automatically satisfied by introducing a 
stream function % defined by 

_ Me Tepe a, soem 
= pe and a cand 

For convenience a new variable R is introduced in place of r; R is defined 
by 


r 
R? = 2 { p*rdr, 
a 


where p* = p/p, = 14+-fe. (8 = m,/m,—1, where m, and mz, are the 
molecular weights of gases 1 and 2 respectively and p, is the density of 
the atmosphere.) 

In terms of the variables (x, R) equation (1) is 


lab ap. ob Oh abory F @ (1 ap 
ee at oh andy —— —.. = aR i d(c)R we ome A He 
ROR ox ER@Rex ox dR? ~ Fal one a) 
In this equation d(c) = p*p*r?/ R® and a = p, when the flow is laminar, 
while for turbulent flow ¢(c) is (p*)*r?/ R? and a stands for py€9. (u* = p/p, 
ft, being the coefficient of viscosity of the atmosphere.) 
Next the equation of diffusion is 
Oy Gc hh Gc é r? 4 


<2 <> > coum — — *)\2 
dRda dz dk") oR 


(4) 


(5) 


*__R 
eR eR 
where o stands for the Schmidt number D,,p,/u. or E according as the 
flow is laminar or turbulent. 


+ The mean value is to be taken in the case of turbulent flow. 
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3. Boundary conditions 
As the boundary of the jet is approached, i.e. as r > «©: u and c tend 
to zero; while on the axis of the jet, ie. r = 0: @u/ér = éc ér = 0. In terms 
of the new variables these conditions are 
L ab 
— —.and ¢ tend to zero as R > x 
ReR 


and w& and pas = 0 when R = 0 and lisa AR a4 = 0. 
4. The first approximations to solutions of (4) and (5) respectively 


The stream function for the round jett 


m u| *) (6) 
I+€ 
which was found by Schlichting (3), is a good approximation when the 
density of the fluid is sensibly constant, that is, at large distances from the 
orifice. 
This approximate form of 4 may now be used in equation (5), to give 
the following equation for the first approximation to c: 


(7) 


eb be eh be é (n) 


= —e ae ome oe UO — 
@Rex cxéR eR\ eR 
neglecting terms of order higher than the first in c. 

The appropriate solution of this equation which is valid at large distances 
from the orifice, is of the form 


c x ma(€). 
The index m is fixed by a condition, which arises out of (5) and the 
boundary conditions, namely, that the flux of molecules of gas | crossing 
a plane normal to the axis of the jet is a constant which is independent 
of x. This leads to the value unity for m. When the above form of c, 
with m = 1, is substituted in (7) the ordinary differential equation 
g i‘ € dg 1 lot eg a 
(i+é)?' 1+€dé * * dé\ dé 

is obtained. The solution of this equation which satisfies the boundary 
conditions g(0o) = 0 and dg/dé = 0 at € = Ois 


g(£) = b(1+€)-**. (8) 


+ The variable é is equal to a*R*/z*, a being a constant whose value depends on the total 
momentum flux from the orifice. 
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The constant 6 in this result depends on the initial flux of molecules of 
gas 1 out of the orifice. It is given by 


M, = 2(1+-B)obp,/(o+2), 
where WM, is the total mass flux of the gas 1 out of the orifice. 


5. The second approximation to the stream function 

Since the first approximation to the variation in density depends mainly 
on the axial concentration b/x, it would appear that the most natural form 
in which to express the stream function is the series 


$= ax “ott F(€)+.. | (9) 


This form of expansion turns out to be incorrect as the following analysis 
shows. 

The expansion (9) is substituted in equation (4) and the coefficients of 
b/x equated to zero. (The terms in ¢ are replaced by the approximate 
value g/z.) This gives for F,(é) the equation 

143, 6 
ot Tee ae 
where primes denote derivatives w a respect to z F, stands for 4¢/(1 
and ¥ is a function of € defined by 


d(c) = 1424+... when ~ is large. 


ha > Zl¥EFo) (10) 


Write P =dF,/dé and change the independent variable from €£ to 
t = (1+ €)-*. Then (10) becomes 
(1 om <-+6P = $, (11) 
where = (d/dt)| V&F5]. 
Equation (11) is of Ss type and its complementary function 
is a linear combination of the integrals P, and P, defined by 
P, = ®(3—4t) and P, = P, log[t/(1—t)|—4—t+ Ye. (12) 


Its complete integral is 


P=-—P, Snare f § 5 P. dt, (13) 


where the lower limits of isteqration may be arbitrarily chosen in the 
general case. Now one of the boundary conditions, corresponding to 
u = 0 on the edge of the jet, is P= 0 on t= 0. This condition fixes 
the lower limit of the second integral in (13) at the value zero. The other 





THE MIXING OF A JET OF GAS. II 397 


condition is that the velocity on the axis of the jet is finite; this implies 

a finite value of P on t = 1. However, it turns out that this condition 

cannot be satisfied since P, has a logarithmic singularity at t = 1 and its 
1 

multiplying factor at t = 1, namely { (S/f)P, dt, does not vanish. 


0 
Thus equation (10) does not yield a satisfactory solution. The flaw in 


the method is that the wrong form of expansion for the stream function 
has been used. Recently some problems in which this type of difficulty 
occurs have been met and solved, notably by Kuo and Ordway (4) and 
by Stewartson (5). Kuo used a modification of Lighthill’s (6) technique 
for rendering approximate solutions of differential equations uniformly 
valid. 

The method, in the present instance, is to expand the stream function 
in the series 4¢ 


b 
<== X —— - F. “T cee ; 9a 

op ot lisetx i(€) | (9a) 
a series which preserves the form of the previous series (9), but in which 
x has been replaced by X and é now stands for a? R?/ X? instead of a? R?/x* 
as earlier. X is defined by the series 


r2=X+ e%(X)+.. 


Since the first approximation to the concentration is now g(£)/X, the con- 
centration on the axis is approximately b/X, and it will be seen that the 
series for y and for x will be useful only when the axial concentration is small. 

The crux of Lighthill’s and Kuo’s method lies in the determination of 
x,(X). The value of 2,(X) is chosen so as to nullify the effect of those 
terms in the differential equation which give rise to the singularity in the 
solution for Fy(€). It will be observed that additional terms arise in the 
differential equation when the derivative with respect to x is changed to 
a derivative with respect to X. Now 


(‘ | ee 1+3("-3 4 1 @ 
dzln dz 0X yao Eaten ae) 


b 2,\] @ 
= l am aa Je ee = 4 a ¥ is 
x(n x)| ax when X is large 


Thus equation (4) becomes 


i eb eh eb Op ap ad| ee 


ReReX ' @R@ROX eX oR? 


=aR ba 
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when the independent variables are changed from (x, R) to (X, R). The 
change of variable thus introduces new terms into the equation. These 
terms are of order (x,—2,/X)b/X times as large as the lowest-order terms 
in (14), while the terms which give rise to the singularity in F(é) are of 
order 6/X times as !arge as the terms of lowest order in (14). The function 
z,(X) is to be chosen so as to nullify the effect of the terms which give 
rise to the singularity. A necessary condition for this is that (2,—2x,/X)b/X 
should be of the same order of magnitude as b/X, i.e. that 2,—2,/X = A 
where A is a constant. This leads to x, = AX log X+ BX where B is 
another constant. The transformation of the coordinates is thus defined by 
- b b 
r= X(14;AlogX+ 5B... 
The equation for F\(€) may now be found by substituting the expansion 
(9a) in equation (14) and then equating the coefficients of b/X to zero. 
|The concentration c is replaced by its first approximation g(£)/X.] This 
equation in F(é) = P has the same left-hand side as (10); but its right- 
hand side contains the term 
ae as wellas — olveri |. 

Thus the equation for P is, using t = (1+-£)-! as the independent variable, 
d? dP 

me; 


where Ss, = pt —2)+8. 


(1—t)— +6P =S,, (15) 


The value of A is to be chosen so that the singularity in the solution of 
equation (15) is eliminated. Since the boundary conditions imposed on 
the stream function do not depend on the x-coordinate they are unaltered 
by the change in variable from xz to X. They are therefore the conditions 
stated earlier. The complete integral of (15) is 

t t 


p= —P,[>P,dt+P, [ 2Pat, (16) 


the lower limits of the integrals being in the general case arbitrary. 
Now the analysis already given for equation (13) shows that the lower 

limit of the second integral in (16) is zero, and that 

Ss 
2 P, dt = 0, 

0 
a condition which fixes A. In the special case when the Schmidt number 
is unity A has the value 12(}8—}(n+8)+ 4n). Taking the still more 
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specialized case of turbulent flow in which n = 28, or laminar flow in 
which 8 = y (which gives the same value of n), A = —}8. When this 
value of n is used 


F4(€) = #P[4t2(3—4t)log t+ 6142— 10004] +a, #(3—4t). 
The constant a, is undetermined because the part /, = ((3—4t) of the 


complementary function of (15) satisfies both the boundary conditions. 
The above expression for F(€) may be integrated to give 


F, = 26[ —4(3t—2#)log t-+-52—4.9¢— 41? +1093] +a, (1 — 3t-+ 2°). 


The stream function of the jet is therefore 


+x hi+| 


Using the inverse relation 


1b b 
X = 2z/)1+-—-Blogz+ B-+...], 
' #8 3 S** zs 
the stream function in terms of the old coordinates becomes 


L x2 4 #8(1—3t+ 21) log? + 
x x 


b 
+ > ise — 4(3t— 2#?)log t+-52 — 49¢— 40? + 10073} 4 


+-a4(1—3t-+-2#?)]4 o(2) (lee) |]. (19) 


where a, is an arbitrary constant. 


6. Conclusion 


The approximate stream function has been found for an axially sym- 
metric jet of one gas which issues into an atmosphere made up of a different 
gas. This stream function is valid at large distances from the orifice where 
the density changes are small. Since the stream function is not analytic 
about x = 00, it was necessary to use Lighthill’s technique, as modified by 
Kuo, in solving this problem. The stream function, however, contains an 
arbitrary constant which is not fixed by the boundary conditions. This 
kind of indeterminacy arises in certain other cases of approximate solu- 
tions of the boundary layer equations when Kuo’s method is used. For 
example, Kuo found an undetermined constant in his solution for the 
flow of a compressible fluid over a flat plate. He suggested that the 
constant might, in this case, be fixed when the full Navier-Stokes equa- 
tions were considered. This is not the case in the problem considered here 
for it is shown in the Appendix that the arbitrary constant is still present 
in the approximate solution of the full Navier-Stokes equations. Thus 
conditions immediately outside the boundary layer do not fix the constant 
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here. Now Stewartson (5) has also found approximate solutions of the 
boundary layer equations in which there is an undetermined constant. 
He has shown, in the case of the wake past a flat plate, that the in- 
determinacy may be removed by taking the flow at the source of the 
motion into consideration. It seems reasonable, in view of the fact that 
the constant is not determined by conditions outside the boundary layer, 
to take the view that Stewartson’s explanation holds in the case of the 
axially symmetric jet. 

Indeed, these arbitrary constants must inevitably occur in the series 
solution, about the point z = 00, of the jet flow. This is because as more 
terms of the series become known so also does the flow of the jet at 
distances nearer and nearer to the orifice become known. As the orifice 
is approached the specific properties of the jet as it leaves the orifice 
become of increasing importance. Now the first term of the series, in the 
present case, only takes into account the most general properties of the 
jet at the orifice, namely, the momentum flux and the mass flux. As the 
orifice is approached from xz = o a point must be reached at which the 
other properties of the jet at the orifice begin to have a significant effect 
on the motion. These additional properties of the jet enter into the solu- 
tion through the undetermined constants. 

The existence of the particular constant, a,, may be inferred from the 
form of the boundary layer equations. These equations are linear in their 
derivatives with respect tox. Thus once a solution, % = f(x, r) say, of these 
equations has been found it is possible to deduce an infinity of solutions 
ys = f(z+C,r), where C is an arbitrary constant, which satisfy the same 
boundary conditions. Furthermore each one of this infinity of solutions 
is valid because no direct reference has been made to conditions near the 
orifice. 

Now the original stream function is of the form 

= ax[R(6)+ Slog("\H(€) +> Hy) +...) 
oe 2 
where = a*R*/x*. Thus the set of solutions which may be derived from 
this solution is 


b 


t b 
b = a(x +0)] Rié,)+ zoel =r) lev + =a H,(€,)+ ~| 


where £, = (a*R*/2x*),_,..¢. When ~ is large this stream function may be 
expanded in powers of C/x to give 


pink na F(6)+ Plog( 2) H4(6)+ 2(H(€) + 5 Me)) +] 
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The function A(€) is in fact equal to 


rhe : 
lim = ehel6)} 
which is found to be a multiple of 1 —3t+-2¢. 

This shows that the most general stream function contains an arbitrary 
constant when only the conditions at the edge of the jet are taken into 
account. 

Since the expansion (19) is not valid near the orifice, it is not possible to 
find the value of a, in terms of the properties of the jet at the orifice 
directly. The constant a, could be found by matching (19) with a stream 
function which was valid near the orifice. Unfortunately such a stream 
function has yet to be found. 
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APPENDIX 


Squire (7) has found an exact solution of the Navier-Stokes equations for the 
flow resulting from a point source of momentum placed in a medium of incom- 
pressible fluid. When spherical polar coordinates (r,0,¢) are used, with origin at 
the source and @ = 0 along the main forward direction of the motion, the stream 
function is 2.5 7(1 — 22) 

eos 
where » = cos@ and d is a positive dimensionless constant. 

When d is small compared with unity this flow tends, inside the core 

1 « p< 1+o(d), 
to that given by Schlichting for the round incompressible jet. Squire’s solution gives, 
therefore, the continuation of Schlichting’s solution outside the mixing region of 
the jet. The possibility thus occurs of making use of Squire’s solution as a basis 
for the study of the flow when the fluid from the source is different from that of 
the surrounding medium. 

Thus the stream function for this case is expanded in the series 


b b) b 
ae} r| fol) + ~ log(—) ful) + Sulu) - (20) 


where 5/r is the concentration on the axis of the jet. This form of expansion is 
fixed by the condition that it must tend to the stream function (19) inside the 
mixing region. 
The equations for f, and f, are found by substituting (20) into the Navier-Stokes 
5002 .56 pd 


= petfo(u), 
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equations of a compressible fluid and then equating coefficients of (b/r) log (b/r) and 
b/r to zero, respectively. The equation for f, is 


(1—p*) fT — (fot 2p) fi +(6+ 3fo)fi = 9, (21) 
where primes from now onwards denote derivatives taken with respect to p. The 


general solution of this equation in f; is a linear combination of two solutions Q, 
and Q,, where (1—ad)* 





Q, = 1430-9) 


Bai Q; 1+p 
5 cee A te 
(a, —p)? 6 (i= )+ 


(a, —p)* a 04 (a —p)* 


8— 20a} + 18a} up — 6a? p* 
(a, —s)® 





and Q: = 
and a, = d+1, 

Now the boundary conditions are that f, and f, should be finite in the range 
—1 <p <1. Thus the required solution of (21) is that f; is to be a constant 
multiple of Q,; Q, is rejected on account of its logarithmic singularities at p = +1. 
The constant multiple in this solution is to be determined by the condition that f; 
be finite in the interval —1 < p < 1. 

The equation for f% is in fact 

(1—p*) fo —( fot 2p) ft +(6+3fo) fo = Ri 
where R, is a known function of f, and fy. The solution of this equation is 
f, = CQ,+ a particular integral of the equation.t The constant C is not determined 
by the boundary conditions because Q, is bounded in the required interval. Neither 
may C be eliminated by consideration of orders of magnitude, becaus. (), behaves 
in a similar way to fj both inside and outside the mixing region. 

Thus even when the complete Navier-Stokes equations are taken into account 
there still remains the undetermined constant. This analysis therefore shows, by 
a process of elimination, that the value at the constant can only be found by con- 
sidering the properties of the flow nearer its source. 


> 


+ Particular integrals of the equation in fj are not in general bounded in the region 
—1< » < 1. A bounded particular integral exists only for a certain value of the constant 
which multiplies ff. 
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SUMMARY 

The effects of uniform high suction on the steady flow of a non-Newtonian liquid 
due to a rotating disk is considered. Series solutions for the velocity components in 
cylindrical polar coordinates have been obtained in descending powers of the suction 
parameter. The solutions obtained are valid for small values of dimensionless 
constant K (= p,w/), where p and p, are the coefficients of shear and cross-viscosity 
and w the angular velocity of the disk. It is found that beyond a certain value of 
the suction velocity, the radial flow of a given non-Newtonian liquid is reversed 
and the axial velocity at infinity is less than the suction at the disk. 


1. Introduction 
THE present work is inspired by the work of Stuart (1) on the effects of 
uniform suction on the steady flow due to a rotating disk. The motion 


of a non-Newtonian liquid produced by a uniformly rotating lamina with 
zero suction has already been discussed by the writer [Mithal (2)] and 


also by Srivastava (3). The coefficients of shear » and cross viscosity 


uw, have been assumed to be constants. As in the case of zero suction 
the pressure has been assumed to be constant for large values of z, r 
remaining finite. Assuming series solutions for velocity components 
in descending powers of the suction parameter a, it has been possible 
to obtain solutions in ascending powers of K (= yu,w/u) where w is the 
angular velocity of the disk. The solutions obtained are valid for small 
values of K and large a. The results thus obtained reveal that above 
a certain value of the suction velocity the radial component of velocity 
becomes negative, which is in contrast to the velocity pattern for a 
Newtonian liquid. The graphs for the radial component of velocity have 
been given for K = 0-1, 0-2 and a = 2, 3, 4, and 5. The paper also con- 
tains tables for the displacement thickness, momentum thickness, their 
ratio and the axial velocity at infinity. 


2. Formulation of the problem 
Let the disk lie in the plane z = 0 and rotate with a constant angular 


(Quart. Journ. Mech. and Applied Math., Vol. XIV, Pt. 4, 1961] 
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velocity w about the line r = 0. The space z > 0 is occupied by a homo- 
geneous, incompressible non-Newtonian liquid. 
The equations of motion are 


ou, 


ou; a ot 
o( cm) pli 


at ex,” 
where p is the density of the liquid, F, is the extraneous force, and t;, the 
stress matrix. 
The non-Newtonian liquid considered in the present paper is that for 
which the stress tensor ¢,, is related to the rate of strain tensor 
Cx = (att in cartesians| 
2 OX; OL 
by [Braun and Reiner (4)]: 
ti, = Wey, +2, 4; j;~,—Pd iy | 
€ 5; = 0) 5 
with the usual summation convention, p denoting the pressure. 
Let the velocity components v,, vg, and v, in the directions of increasing 
cylindrical polar coordinates r, 8, and z be represented by the relations 
v, = rf (2) 
Vg = rg(z) (2.3) 
v, = A(z) 
and as in the zero suction case (2) the pressure distribution is assumed 


given by p = x(z)+r°d(z). (2.4) 


With these values of the velocity components and pressure, the equa- 
tions of motion (2.1) under the transformations z = (v/w)*¢, f(z) = wF(f), 
g(z) = wG(C), A(z) = (vw)tH(f), and x(z) = —pywX(f) become 


F?—G?°+HF’ = F"—}K(F"?+3G"+2FF’), (2.% 
2FG+HGQ’ = G"+K(F'G'— FG’), (2.4 
HH’ = H"+3KH’'H"+X’, (2.7) 
H’'+2F = 0, (2.8) 


where v = p/p, v, = p,/p, K = p.w/p and primes denote differentiation 
with respect to ¢. 

Since H is assumed to be independent of r, the equations of motion 
(2.5), (2.6), (2.7) and the equation of continuity (2.8) are applicable only 
to the case of uniform suction through the disk (say w,), H taking a 
constant non-zero value at [ = 0. 





ON THE EFFECTS OF UNIFORM HIGH SUCTION 
The boundary conditions are 
F(0) = 0, F(a) = 0 ) 
G(0) = 1, G(a) = 0 . (2.9) 
H(0;} = —a, H(w) =a Tee 
We now define a suction parameter w, R'/U,, where wz, is the velocity 
of suction, R a Reynolds number, and U, a representative velocity. In 


our case, since wy = a(vw)t, U, = rw, R = r*w/v, a is the appropriate 
suction parameter. 


3. Solution of the problem 

In this section we proceed to obtain series solutions in inverse powers 
of the suction parameter a. 

We consider the equations of motion (2.5), (2.6), and the equation of 
continuity (2.8) under the boundary conditions (2.9). On physical grounds 
H may be assumed to be nearly constant for large values of a, and conse- 
quently H’ = F = 0. Following Stuart (1) we get as a first approximation 

G = e-%,. (3.1) 

This suggests that a solution in descending powers of a may be obtained 
and that af is more suitable than ¢ as an independent variable. We there- 
fore write 


n= al 


and transform equations of motion (2.5), (2.6), and the equation of con- 
tinuity (2.8) accordingly. We then have 


F?_G?4+aHF’ = a?F"—ja®K(F?2+3G'242F F’), (3.2) 
2FG+aH@’ = a®@" +a?K( F'G' — FG’), (3.3) 
aH'+-2F = 0, (3.4) 


where primes denote differentiation with respect to y, and F, G, and H 
are now functions of ». 

As in the case of the Newtonian liquid considered by Stuart (1), we 
assume the following solution 


F = Sa~F{n) 


x 
r=0 
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When a is large and K small, series expansions for 7, G, and H may 
be found by the method of Stuart (1) with the following results: 


PF = 2K (e-*7—e-1) 4 sea (€"—€-") + BK 20-4 70-99 + ie “27 $e-n) + 


"2 
+ Efe 4e-24 (94 9n)e- 4 tne ues 


bak (he Hte-#9-+ He e-9+ Ife en) 4 


K 7 2 
+ aaltie Fe + (e—Inde-+ (H+ Bn)e-7} 4 


K* ‘ , 
{—Bille-"-+ gen (9998 tn )e—tn + (8991604 )e 27 


* Tea : 


— (Si — lin)e 274 (15 ay — Sn )e “m\ 


l eT ke , 
+ qt—ie A+ fe-81-+ (G+ ane“ — ($8 4+ dye} + 


"6 "3 "2 3 
+0(K°, = , a , x a) (3.6) 
a*° a a a 


, 


“ ee 
G = €14 bKMe-"—e-M) + — fhe 37 (8 3n)e-) 4. 
+ BKM —be-P$ fet —to-904 e-) + 
l s 1 3n l n\ 
+a — qe" +- (5 —$)e-7} + 


kK3 


+ 3902 


{118¢ “$n 13e-47 + (47 - 9y)e-*9—(8—in)e m4 


5 


+ sigK*(—188e —77 +-194e-67 __ 1985 51+ 30e-41 — 88-37 j he “) 4. 


Kk? . 
5s : j ‘ . 59 \o—m_ 
+ gga tebe "+ Be + (ty +- Bn e-97 — (2552 — 91.) — 99 2m) 


K £ ” 3 9 
n 3 o— $7 »- : 3 343 2)p—m 
7 tino"  — ge *— ($8 -+-4)e "+ (188 — 83 — 2 y*)e G+ 


a rs ' eat 
tag —slae "+ dee +-(h+ de 7+ (— + Hn + dn t)e-1} + 


ve K6 K4 K3 k2 
+O( RY es Ta ae a) 
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and 


a 


me Spe Pe 1 3 ¥ 
H = ore 5 te "—¢ "+4)+4(—te *7-+-e-"—4$)+ 


+i at (4e-* —je-9" +-18¢-29 — Je-1 4-4.) + 

a 

a) —4n 1 379-3 23 -2 7 ~7 1.178) 
a | tje-*9 + Fe-*9 — (B—On)e-9— (H+ One" ++ 


9 K® 
+55 — (the —Be-*+ et e-P9-+ He he-— J) + 


K 
4+. — £13 9—4n__ 49-39 (29.4 ~2n_} __ 3010 
a®' (sue * — fe — (38 +- $y)e-* + (H+- 9 )e-” — Bh} 


K* . —_ 
ae T= 3{— — 938%¢ ~6n 4 19M e- 57) — (3980 io Bn etn + (2587 Se $y )e-o9 — 
— (%'— 15n)e 21+. ($278 —13y)e—7 — Sees | + 


] P 
de ~ { — ie" ak he -3n +(}3 _ $n)e-27— ($$ +- kn)e ~7 + 30h} ol 


"7 6 "3 €) 4 
a=, =, = 2 a (3.8) 
acees € 
It may be remarked here that the values of F, G, and H given in (3.6), 
(3.7), and (3.8) are valid only for high suction, i.e. when a > 2. The case 
when a < 2 will be considered in a later paper. 


4. Some standard parameters 

In the case of three-dimensional boundary layer flows it is preferred to 
define two displacement thicknesses, in two perpendicular directions. In 
our case these two perpendicular directions may be chosen as those along 
the radius and the tangent to the radius. In the present case, however, 
there is no radial flow either at the disk or at infinity; consequently the 
radial displacement thickness has no meaning. The tangential component 


of flow gives - 


5* = | ("2) dz, 
J \rw 
0 


which under the previously defined transformation becomes 


8* = (v/w)t | @ dé. 


Similarly, we have the momentum thickness 


= (v/w)? { aa—@) dt. 
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The ratio of 5* and @ may also be defined by 
H* = 8*/0 (4.3) 


(* is introduced to avoid confusion with the velocity function H). 

For a > 2, we have, after substituting the value of G from (3.7) and 
integrating, 

] 1 MA. .4 ee ee 

Reh an 1 R84 a. KA+..,}, 4.4 

ee | +79" 1 78 at 9a 3200 ie 

Of 7 51 », ; 


O(w/v)* = w IES he, a 4.5 
lop = sa +3 72 a? 36a* 9600 (4.5) 


and the ratio of these two gives 


1 5 K l a 
ae a a ee Big eeeee K*+... 1. 4.6 
' +94% 72 a’ 36a 3600 — 
There is some error in Stuart’s calculations. The values of @(w/v)' and 
H* given by him are 
O(w/v)t = “A! - +, 


2a\ sd 2a 
and H® = (1+ ‘ 
36a* 


We find that these values for a Newtonian liquid should be 


Ww/v)t = “(1 Ad + (4.7) 


~ 36a! 


2a 


and Hf? = (1+ sat) (4.8) 


The following two tables give the values of 5*(w/v)', @(w/v)#, and H* 
for different values of K and a. 


TaBLe 1 (K = 0-1) 





3 4 
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Also, the axial velocity at infinity (say H,,) is given by 


‘ 3 391K , 20 
BE. sen eth n Cis tintin 391 K 01 


ro AERA cnn A) PR ED 4.9 
4a 2a* 64a 288a5 88a? * On 


It can easily be seen from (4.9) that under our assumptions about K 
and a, viz. K < 1 and a > 2, H,, is always negative. Further, when K 
is small and a sufficiently large, we find that the axial velocity at infinity 
will be less than the suction at the disk provided 


Ka* > §, 


(4.10) 
y 
Thus for a given non-Newtonian liquid if the suction velocity is large 
enough to satisfy the above condition the axial velocity at infinity will be 
less than the suction at the disk. The following table gives the value of 
the axial velocity at infinity for different values of K and a. The values 
at K = 0 are included for the completeness of the table. 


TaBLE 3 (—H,) 





K 


orl 








2°062 
2°999 
3°989 
4°989 

x 














5. Discussion 
If we consider the expression (3.6) for F, when XK is sufficiently small and 
a sufficiently large, such that only the first two terms in (3.6) dominate, 
we find that F becomes negative for 
Ka® > § 
wiv, 2 
— —_— Ps a 
v2 3 
From this we conclude that for a given non-Newtonian liquid, reverse 
radial flow can take place for a high suction velocity. Comparing this 
with (4.10) we find that reverse radial flow takes place if the axial velocity 
at infinity is less than the suction velocity at the disk. If, however, 
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the disk acts as a centrifugal fan as in the Newtonian liquid case con- 
sidered by Stuart (1). These facts are clearly brought out from Figs. 1 
and 2 which have been drawn at magnification 100. 


-00«F for K=-1 











—4100 x F’ for K =2 








The velocity component @ behaves in a manner similar to that of the 
Newtonian liquid. 
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A PROBLEM ON THE SLOW MOTION OF A VISCOUS 
FLUID BETWEEN TWO FIXED CYLINDERS 
By K. B. RANGER 
(Bedford College, University of London)t 
[Received 10 October 1960] 


SUMMARY 

The two-dimensional Stokes flow between two fixed coaxial cylinders due to a 
line source and a line sink on the outer boundary is determined in terms of an 
Earnshaw stream function. 

Some numerical values are given to illustrate the variation of the force per unit 
length on the inner cylinder with its radius, the outer boundary of radius unity 
being kept fixed. These results are compared with corresponding values obtained 
from the irrotational non-viscous problem. It is found that in the transition from 
a very viscous fluid in slow motion to a non-viscous fluid in irrotational motion the 
direction of the force is turned through a right angle. 


1. Introduction 
RAYLEIGH (1) has considered the two-dimensional Stokes flow within a 
circular cylinder due to a line source and sink of equal strength situated 


at opposite ends of a diameter. The fluid 
enters and leaves the enclosure by means 
of two small perforations made at oppo- 
site ends of a diameter, the boundary 
being impenetrable elsewhere. 

The present work extends Rayleigh’s 
problem by considering the flow between 
two fixed coaxial cylinders due to a line 
source and sink located at arbitrary 
fixed points on the outer boundary 
(Fig. 1). It is noted that problems of 
flow between two cylinders are not 
in general unique. This is due to the 
possibility of obtaining solutions of the two-dimensional biharmonic 
equation corresponding to a circulation in the region and which also 
satisfy the condition of no slip on the boundaries. Uniqueness in the 
present problem is achieved by taking the circulation in the flow to be 
zero. The motion is prescribed in terms of the Earnshaw stream func- 
tion and found in the form of an infinite series. 


+ Now at Department of Mathematics, University of Toronto, Toronto 5, Ontario, 
Canada. 


Fie. | 
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The main interest is to determine the force per unit length exerted by 
the fluid on the inner cylinder. Considering a cross-section of the motion, 
the line perpendicular to the source and sink through the centre is a 
line of anti-symmetry for the velocity field. It follows that the stresses 
round the inner boundary are themselves anti-symmetrical and the 
total force acts perpendicular to this line. Calculations show that the 
magnitude of the force per unit length of cylinder is proportional to 
the sine of half the angle subtended at the centre by the chord joining the 
source and sink. Numerical values show that by fixing the radius of the 
outer boundary at unity, the variation of the force with the radius of 
the inner boundary increases rapidly when the latter exceeds one-half. 

The flow is compared with the corresponding irrotational non-viscous 
problem. Here the pressure, being proportional to the square of the fluid 
speed, is symmetrical about the line of anti-symmetry for the velocity 
field. Further the maximum fluid velocity in the region occurs on the 
mid-point of the line joining the source and sink, whilst the minimum 
fluid velocity occurs on the line of anti-symmetry and on the side of the 
inner boundary remote from the source and sink. Thus the total force in 
this case is directed parallel to the line of anti-symmetry towards the 
mid-point of the line joining the source and sink. The couple on the inner 
boundary, is zero as is the case in the viscous problem. 

A correction due to Basset is applied to the pressure, which in effect 
partly takes account of the neglected inertia terms from the equations 
of motion. However, the results are not affected by this correction and 
it may be concluded that if the motion were slow enough or the fluid 
viscous enough, the force component perpendicular to the line joining the 
source and sink would vanish. 


2. The equations of the problem 

It is first convenient to determine the Earnshaw stream function for 
the Stokes flow within a fixed circular boundary of radius a, due to a line 
source and sink of equal strength m situated at arbitrary points on the 
boundary. 

Considering a cross-section of the motion, the centre of the circle is 
taken as the origin O of polar coordinates (r,@), and the coordinates of 
the source and sink are A(a,«) and B(a, —«) respectively. 

The stream function for the combination of source and sink alone is 

pga m| tan oe acon) Poo §+asin *)} 


- — . (1) 
r cos @—a cos « rcos 0@—a cos a | 


This is a streamline for r = a, except for the points of discontinuity 
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@ = +a where the flow corresponds to a source and sink; and for r < a 
may be expanded in the form 


a 
3 - 1/r\" . 
jy = —2maxa—2m > —|—} cos n@sin na. 
1 n\a 
na= 


The complete flow may be represented by 
$= fot, 
where y, satisfies the equation 
Vi #, = 0 
ge 2 o2 
wi m=(Z+12+42) 
In the present case, an appropriate general solution of (4) is given by 


%, = s (A, r®+ B, r®+*)cos n@, (5) 
n=1 ' 


A,, and B,, being arbitrary coefficients. 


3. Boundary conditions 
The components of fluid velocity (u,v) referred to polar coordinates 
(r,@) are given by 1 ay aus 
u= —-—, =z: 
r 00 er 
u and v vanish over the boundary r = a except at the source A and sink 
B. Thus the condition to be satisfied by the normal component wu is 
¥,=0 at r=a 
and by the tangential component », 
_ ae 


—¥. at r=a. (7) 
or 


Substituting (2) and (5) in (6) and (7) we obtain the following system of 
linear equations to determine the coefficients A, and B,,. 


A,,a"+B,a"+? = 0, 
nA, a"+-(n+2)B,a"*? = 2msinna 


Mm sin na 
“an 


giving A,=-— and B, = 


Then ¥, is given by 


) > (Z)"sin na cos n6, 
a 


n=1 
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and #% can be written as 


= —2ma+ > Ja(roos nd (10) 


n=1 


where f,(r) = —2m m(2)'{5+5(1-3 }sin na. (11) 


\n 


Again by differentiating (2) with respect to r, we find that (9) is equivalent 
= ( - alire 


to 


and hence 


t= bo+4(1—S)r rte 


or 


Aa ES 


+4mar(1— 


rsin @—asina _,/rsin @+asin « 
— tan-}{ —__—— 
r cos 8—a cos « r cos 6@—a cos a 

\ (r?+-a?—2racos(0+-a))  (r?+-a®— 2ra cos(0—«))|" 
It is readily seen from this expression for 4, that the flow around A and 
B corresponds to a source and sink, the fluid velocity vanishing at all 
other points of the boundary. 








‘) sin(@+-«) sin(@—«) | 


4. Flow between two coaxial cylinders 
A cylinder of circular cross-section and radius b (< a) is now introduced 
into the flow and is fixed in a coaxial position with respect to the outer 
boundary. 
The stream function for the motion is written as 
b= fotbittes 


where yf, is an appropriate general solution of (4) given by 


=) nl (r)cos n8, 


n= 


where 
g,(r) = A,r"+- =2+0, usa De (n > 1) 


g(r) = Ayr+ByP+ [4 Dyrlogr, 


4, B,, C,, D, (n = 0,1,...) being arbitrary coefficients. 


5. Boundary conditions 

The conditions of flow within the outer boundary are satisfied by the 
stream function %,+¥~,, hence it is sufficient that the velocity components 
derived from ys, vanish at r = a, i.e. 


8 treme. (15) 
cr 
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On the inner boundary the fluid velocity vanishes everywhere, i.e. 
op 
= —_- 0 at T= b 
or 


so that using (10) and (13), (15) and (16) are equivalent to 


9,(4) = g,(4) = 0) 
fx(b)+9n(b) = fr(b)+-9n(b) = 05 


where the dash indicates differentiation with respect to r. 

Equations (17) constitute sets of simultaneous linear equations to deter- 
mine the constants A,, B,, C,, and D,. These are as follows: 
D, 


n - 0 
a” 2 


(17) 


b 
n 


A,,a" + Bn o qnity 
a 


nB,, Y —2 
F +2 
nA, tees + (n+ 2)C, a®+? — ~— d 


ee D, 
A, "+5240, b+ <8 = —f,(b) 


—2)D 
nA, 2 oe 4+ (n+2)C, ones (0S) 


and for n = 


Aya+Bya4© + D,aloga = 0, 


Y 


( 
A,a+3B,a*—— + D,a(loga+1) = 0, 
1 1 1 £ 


A,b+B, 44 D,blogb = —f,(b), 


A,b+3 Bp + D,b(logb+1) = —bf1(b). 
bh 
Solving we find 


b\" . 
= —mb"-*{_]} sin na x 
a 


€ 2 
x (HS) 4 2)b?(a2" — b2”) - n*a*"(a?—b?)}+- 


n 


oii +9 b? 2,,2n/ 72 h2 __ h2(q2r __ p2n | 
-(n+2))1— n®a2"(a? — b?) —b*(a 52") 
a® t yy 
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2k, B, = mbr-+()"sinnax 


<(t! 


. =a) (at —beyn-+-m(n—2)(a*§—62}0™™-4} + 


+(n+2)(1 ‘a ca) tat at—08)—(a™—o™)| 
2k, C, = maror-+(")"sin na x 


x (Ce - ia) (at—be)(ntat—(n— 1)(m-+-2)b?) —n*a?"(a?—b*)}+- 


+(n+ 2 — 5} natn(at—2)+ (a —b) (0 _ 1)b*—na)} | 


b\" . 
2k, D, = —mar-ron-+) sin Na X 


x Ce a) a2"—b™)n((n-+1)a*—(n-+2)b%)—n(n-+-2)n?(a2—b2)} + 


2 
+(n4.2)(1—}fnat(at—b2) + (a™—b%)(b¢n—(n-+ 1a} 
for n > 1; where 
k,, = {(a2®—b)2— n2(a2—b2)2q2n-2h2n-2) 
For n = 1, A,, B,, C,, and D, are given by 
k, —mMsin a 


1 ga(a— BP) * 


‘ (a( — 5) |(2loga+ I)(at 0) — dab log ?— 20%(a*—b2)} + 


+(3— 5) log a+ 1)(a?-+ 36?)(a?—b2) + 4a*b? log” — 2a%(a*—b*)| , 


b? si ae 
k, B _™ (sal 


= | 
- Misia 2 . 2 om 
1 = Fat) 3 b?) + 2a* log | 


f | 
iuicel sie inci 2__f2 2a2 —)}), 
( (2 +-2a* log 
k,C, Mm Sin a 


b? b 
= 2a(a*—b) (s—3){e log = _ Har—b)| _ 


_ a(t — 5) at—0+ 2a? log 7 


a2 
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where k= (a*+b4)log” + (a*—b). 
For large enough n it can be shown that 
A,r*| = = Cnt?) < fin)(3) 


\B,| _ 0*| D, b\ 

re; Pega Mee n)i — 

| = alata] <tm(2) 

where f is a polynomial in n of degree not greater than 4 with positive 
coefficients. These inequalities suffice to infer that ¥, is uniformly con- 
vergent in the region of flow. Again it is noted that the stream fvaction 
~ for the flow is unique in virtue of a theorem given by Lord Rayleigh (2). 


and 


6. Calculation of the force on the inner cylinder 
The components of force (X,Y) on the inner cylinder are given by 
20 
Xi+Yj = ( (R,),-, 6 d@, (18) 
0 
where i and j are unit vectors directed along the axes 6 = 0 and @ = 4r 
respectively, and R, is the stress vector, which, in the usual notation, is 


iven by 
given by rR, = —pr+pf(r.v)qa—q+V(q.r)}, (19) 


where , is the coefficient of viscosity. In particular the stress vector over 
the boundary reduces to 


R.=!/—pri+pire! , (20) 
. | ér® Jrmb 
where r and 6 are unit vectors directed along and perpendicular to the 
radius vector respectively. 
Substituting (20) in (18) and integrating the pressure term by parts, 


we find o 
a2 
X= [ op _ Y\ bcind dd 
26 Or* } <5 





= f (2 _. a] bcos 6 d8 
06 r=b 


The equations of motion referred to polar coordinates are 


Cp _—# oe ws, lop _ 29 
ori 26 ("1 #), + 26 a EY) 
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giving in particular 


op, _ (Oh LM pe 
@)_.- later. 


Substituting (22) in (21), the final expressions for the components of force 


are given by Pe 
> Ay , 
= 2 — es 
X = pb | (3 ) sino dé 
0 





20 
‘ Y (Ay 
Y= — a | (55) cos 6 dé 
e : or® r=b 


) 


Replacing % by (10) and (13), (23) becomes 


27 


x. >t f2(b)+97(b)} | cosnésin 6 28, 
sot 0 


27 
Y = — uit > { f""(b)+-9%(b)} | cos n6 cos 6 dé 


0 
showing that the evaluation of the force components depends on the inte- 
grals an Qn : 
( cosn@sin@ dé and ( cos n@ cos 6 dé, 
0 0 
the former of which is zero for all positive integers, and the latter being 
zero for n > 2 and equal to zw for n = 1. Thus 


X=0 and Y = —mpb{f7(b)+97(6)} 
which on calculation is 


Y= 4rmy(3a?+-6*)sin « (24) 
a{(a?- +b)log(b/a)+a?— ~b?}" 5 

The denominator of this expression is negative, indicating that the total 
force on the inner cylinder per unit length acts coparallel with the line 
joining the source and sink. 

By letting m > «©, «0, so that lim 2maa = M,, the source and sink 
tend to a doublet of strength M, at (r = a, @ = 0), and directed along the 
tangent to the boundary at this point. 

The force per unit length of cylinder in this case is 


Se. (3a?+- b?) 
a*{(a*? + b*)log(b/a)+a?—b?}" 





Y = 27M, 
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Some numerical values are given below for the variation of the force with 


the radius of the inner boundary, the outer boundary being fixed with 
constant radius unity. 





$ i 








47m sin a 




















mite "7 "35 121°82 | 296°88 





The couple on the inner cylinder is given by the formula 


2a 
G.k = { (rAR,),.»b d0 


0 


22 


which gives G = pb | (=) dé, 
r=b 


er? 


k being the unit vector perpendicular to the plane of motion. 
Using (10) and (13), (26) becomes 


27 
G = pb? 3 (f7(b)+-9%(b)} [ cos nd a0 
n=1 Py 
which is zero. 

Finally it is of interest to consider the change in the force per unit 
length on the cylinder when the flow takes place with an incompressible 
non-viscous fluid. 

The stream function, which is easily determined by the familiar methods 
(3) of potential theory, is given by 


,(7sin @—asin « _,(7sin 0+asin « 
ub m tan~!{ ———_—_—_—_—__ ]} 4 m tan -*+| ——______ 
r cos 0—a cos « r cos @—a cos « 


x 


1 b2 n yp2n —_q2n : 
+2m > -_|—} ———— cos n@sin na. 
— niar} a*”—ph” 


n=1 


The force components are found to be 


X = Y = 0, 





Srm2p Ss (=\"" sin na sin(n+- 1l)a 


b a a (a2” — §2")(q2+D_— prn+D)’ 


p being the density of the fluid. 
In the limiting case when the source and sink tend to a doublet of 
strength M, at (r = a, = 0) and directed along the tangent at this 
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point, the force is given by 


x= BAe S (2 n(n+1) 
si ( 


ba? a a2” — h2)(q%n+1)__ pXn+1)° 
= 





Some numerical values are included below to show the variation of the 


force with the radius of the inner cylinder, the outer = fixed with unit 
radius: 





b 3 i 
10°X 
aMip : "76 | 138157 | 683-39 
=e 





























The force is thus directed towards the centre of the line joining the source 
and sink, increasing rapidly when 6 exceeds one-half. It is easily verified 
that no couple acts on the inner boundary. 

The direction of the force is turned through a right angle in the transi- 
tion from a very viscous fluid in slow motion to a non-viscous fluid in 
irrotational motion. It is noted that the stress distributions round the 
inner boundary in the viscous and non-viscous problems are respectively 
symmetrical about the axes 6 = 47 and 6 = 0. Again the couple on the 
inner boundary is zero in both cases. 

In a modified theory of Stokes flow attributed to Basset (4), the pressure 
term p is replaced by p+4pq*, where q is the fluid speed at a point deter- 
mined from the approximate equations of motion. The theory, in effect, 
is an attempt to partly take account of the neglected inertia terms in some 
regions of the fluid where the Stokes approximation is invalid. In the 
present problem this may be the case in small regions surrounding the 
source and sink on the outer boundary. Applying Basset’s correction, 
the force components (X,Y) per unit length on the inner boundary are 
increased by increments (X’, Y’) “ne where, from (21), 


= x =e fe Cin 6 db 


and —}pb fae ) o08 0 dé. 

0 
However, over the inner boundary @(q?)/2@ is zero, so the resultant force 
is unchanged although the pressure distribution is no longer symmetrical 
about the axis @ = 47. Again from (26), the contribution from the pressure 


vanishes everywhere on the boundary, hence the resultant torque is still 
zero. 
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THE MOTION OF VISCOUS LIQUID PAST A 
PARABOLOID 


By D. J. MATHER 
(Atomic Weapons Research Establishment, Aldermaston, Berks.) 


[Received 27 March 1961] 


SUMMARY 


An approximate solution for the steady flow of incompressible viscous liquid past 
a paraboloid of revolution is described. An assumption is made for the form of the 
Stokes stream function and substituted into the Navier-Stokes equations using 
paraboloidal coordinates. After making suitable approximations, a non-linear 
differential equation for a function F is deduced. The solutions of this equation 
depend on the Reynolds number of the flow considered. Examples found by 
numerical integration are given to illustrate the properties of the function F for 
Reynolds numbers varying from 0-0001 to 10,000. An expression for the pressure 
is found, and it is shown that this approximate solution tends to the perfect fluid 
flow away from the boundary, allowance being made for the displacement effect of 
what may be called the boundary layer. 


1. Introduction 


THE steady flow of viscous incompressible liquid past a paraboloid of 
revolution is considered. The method used by Dean (1) in considering 
the flow of viscous liquid past a parabolic cylinder, suggests a similar 


approach in this case. 


2. Egnations 
The Navier-Stokes equations for the steady motion of a viscous liquid 
are made non-dimensional by writing 


2’ = Wejv, y' = Wylv, 2 = W2)v, 
where 
y2+2'*=-7'8, p= pW’, u= We’, v= We’, w= Ww’, (1) 


and W is the undisturbed velocity in the x-direction. 
The Navier-Stokes equations whose solution is required become, 
dropping the dashes, 


grad{p+4q7} = q A E—curl E, E = curlq 
together with the equation of continuity 
divq = 0. 
Paraboloidal coordinates are defined by 


2’ = &2'—7?, r’ oo 2&n (» a 0). 
(Quart. Journ. Mech. and Applied Math., Vol. XIV, Pt. 4, 1961] 
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The motion is axisymmetric and there is no rotational component of 
velocity. Resolving the velocity into components (U,V,0) in the direc- 
tions (€, n,¢), the equation of continuity (3) becomes 

sg ltlE+ aU} + = Ent 2 V} = 0, 
which may be integrated by introd~’ 1g a stream function y’ defined by 
i <7 oy V= ; se 


~ REEF) On?” ARDEP) CE " 


e fv (€?+ 2 ab (+7)U}| =, say; 


t= se; Hag 
ie. (3 i SAT soit ol ey Ley’ 1 yp’ 


8én(€*+-m)| 2€? “On? € OE y OnJ 
Equations (2) can then be written 
P_Q M2 Q P_ Aw Mo 
o& 2 0 OH 9 én = - Eq On é 
where P= p'+4q”. 
Eliminating P from equations (7) gives 
ef 22, 22), JO_O), fe2_O 
ee" on? "OE ES" leq” af 
lap’ Qa’ Q)__Q ay’, O ay’ 


2| 0 Gn Gy CEJ 2m OE * 2 ey” 
3. Method 


It is now shown that an approximate solution of (9) is given by 
i = —€F(y—a) (a> 0). 
With this expression for %’ the velocity components are 


i "Any (€+ 7)’ 
The reduced vorticity Q is given by 
€  {F'(n—«)_ F"(n—2)| 


(+97)? “Wes 3 
The condition that (10) should be an exact solution is 
2F” = 3F" oo: 2F” red {ie <2 | F’) 
Fiv—_—__4-—__——,_+F - +(\F mz 
.s i> oo wll ha 
F’ 


* ele ral )enF \ 

= | Fv 2p Fn PP" +| P’—=\(2nF’— F—1) 
+x)" . a)" zB 
+ gapgaal PoP M42 F 9 F)— 4 PY 
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We now assume that J satisfies the left-hand side of the last equation 
put equal to zero, i.e. that F satisfies 


pe_2P" SF" 3F' FF" 3FF’3FF FF’ PF? 
ners ee nea ap ete er ore ate pete Sass 
7 I 7 ¥ 7 7 7 7 


This equation can be integrated once to yield 
7? F” —yF" + F'’+7FF’—FF' = 0, (14) 


where the constant of integration has been taken to be zero by virtue 
of the boundary condition on F for large (y—«), discussed later. With 
this choice for F, the error term £ on the right-hand side of (13) becomes 

simply 
= pe satel | 
8(&+ n?)?| 


For given », EZ takes its maximum value £,,(y) when £ = ». This gives 


P(r ——\+ F(R nF"). (15) 
7 


32y*hy = (er ——)o -9F’). (15) 


This error is tabulated as a function of » with the other results. 
The conditions that the velocity should vanish on the surface of the 
paraboloid are, from (11), given by 


F(0) = F’(0) = 0. (16) 


And since ys’ = —2£*»? is the stream function for the undisturbed inviscid 
stream, the condition on F for large (y—«) is 


F’(o) = 4. (17) 

The surface 7 = « is the paraboloid 
r? = 4a(x+a) (18) 
a? = Waly = R. (19) 


It may be noted that, when a = 0, the function F = 27? is a solution 
of (14) which also satisfies (16), (17). However, the velocity component 
U does not vanish on » = 0, so this does not give a solution of interest. 

Equation (14) has been integrated numerically for several values of «. 
Beginning with an assumed value of F”(0), step by step integration was 
carried out, using a Taylor expansion method, until a value of F”(a) 
could be deduced. If this did not yield the correct value, 4, another choice 
of F"(0) was made, and the process repeated until agreement was reached. 
The results of several such integrations are tabulated below, together with 
the asymptotic properties deduced from the numerical work. 

The form for the pressure corresponding to the expression (10) is now 
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found. Substituting in equations (7) the following two expressions are 

found: aed F’| 
o& — -4(€#+- 9%) aS 
ge (rr | (F’ 
én - weTFA 7 er h 


where (14) has been used to obtain (20). The integration of (20) with 
respect to & yields 


(20) 


(21) 


lL {pw FF 
= —.__— | F”— —_}+-G(n), 
8(€*-+ 7)| > Tih 
where G(n) is a function of » only. Equation (21) is satisfied with this 
form for P if 


na AP. PM | eS , »)\ 
gout 5 r|+ aera tra wl 1 


This equation is approximately satisfied by ignoring the second term on 
the right, which may be compared with (15), obtaining 


F” 
G(n) = +¢. 


327? 


This gives Pie ae ie 
_ = 329 t ETM 


sak =| +C. 
i) 


The reduced pressure can now be found from (8) and (11) in the form 


; ee Se Pe ty | : 
SS | menial — — « ——}+C. 23 
8(€+- n*)| ee ee nr 
From the numerical results it is seen that the properties of F(n—.«) 
are, for large (7—«) 


F = 2(7n?—«"?), a . PF" = 4, (24) 
where «’ is a number whose value is close to that of «a. It is now shown 
that %, the actual stream function, given by 


9,2 
= =-¥ = = -= (4? —a')E* 


is the correct expression for the irrotational flow of inviscid liquid past 
the paraboloid, making allowance for the displacement effect of the viscous 
layer near the boundary. 

The equations governing the flow of inviscid liquids are 


curlq = 0, divq = 0, 
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subject to the condition that on the boundary the normal velocity com- 
ponent is zero. That is, the solution of 


Oy a 1 op 1 oy 


dently enim ancients slits pce mata ame 


08 oni §; , om, 


is required, where 
a ow Oe oyp a 1 opp, 
S48; mV (G+ nh) en,’ * 4; ma (G+ nb) 
are the actual velocity components in the directions £,, 7, increasing, 
x+ir = (€,+%,)*; 





V, satisfies the conditions 
y=0 on =” V,>—W for large »;. 


no is the parameter corresponding to the surface of the paraboloid, so that 
from (18) it is deduced that n, = a. It may be verified that the required 
solution is b= —2W Ey? —7}). (25) 

Hence apart from the slight difference of «’ from «, the stream function 
derived from (10) tends in the limit of large (y—«) to the perfect liquid 
flow, this difference being due to the viscous effects at the boundary. 

The pressure of the viscous liquid is pW*p’, where p’ is given by (23). 
For large (»—«), taking «’ as approximately equal to a, the limit is, apart 
from a constant, 

26+) nf 2UG+nDl aid 
and (26) can be derived from (25) and Bernoulli’s equation. 

There appears to be no upper limit to the value of R in the integration 
of equation (14) for F. It may be noted, however, that no solution exists 
for a slow, i.e. Stokes, flow when R = 0. This may be seen as follows. 
If the quadratic terms in (14), arising from the non-linear terms in (2), 
are neglected (which also means E = 0), then the resulting linear equa- 
tion can be integrated to yield 


F = e+fn?+977 log », 


e, f, g being constants. But with this form of solution it is impossible to 
satisfy the boundary conditions (16), (17). 


(26) 
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a = 0-01 R = 0-0001 
F FF’ 

0000 0000 
OOo! 0-048 
OO17 0-298 
0-063 0°635 
0146 1°023 
0-269 "445 
0-436 1-888 
0°647 2°345 
O°go5 2°808 
1*209 3°273 
1°559 3°736 
2-398 4°645 
3°416 5°525 
4°606 6-378 
5°965 7°209 
7°489 8-026 
9175 8834 
11*022 9°638 
13°030 10°439 
15*198 11*240 
17°526 12°040 
20°014 12°840 
22°662 13°640 
25°470 14°440 
28-438 15°240 
31°566 16°040 

The form of F for large 7 is 

F = Xy*—a’¥), 
where a’ = 0-444, 


‘ 


§ 


"000 
0°599 

1°248 

1°939 

2°660 

3°400 

4°147 

4°888 

5612 

310 

6-976 

8-199 

9°279 

652 10°244 

791 II°131 

102 11°973 

57° 12*7990 

*217 13°597 

016 14°399 

*976 15°200 

096 16°000 

376 16°800 

*816 17°600 

416 18-400 

“176 19°200 

096 20°000 

The form of F for large 7 is 

F = %(y*—a), 

where x = 1-566. 


a Hue 


sO DOew se 


PWIA DIM 
won ns 
wv 


fat ott tt ttt 
Saunt wWnw 6 3 


oO 
~ 


. SS oe 
of & 
a] 
an 
Bd 


Qe Sw OwmnewnNdNen 00009 
as “~ 
SES £3 
te & & yO 


ewww nnn ee 


ws ~~ Ww 
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a= 10 R= 100 
F F’ Fr —Ey 
"000 0"000 39°812 0°000 
"200 4°000 40°182 0-488 
o8o1 8-031 40°390 0955 
1°806 12°068 40°286 1°391 
3°214 16°073 39°735 "774 
5018 19°997 38°629 2-081 
7°208 23°779 36-098 2°291 
9°767 27°358 34°573 2*387 
12°671 30°675 31°689 2°365 
15°892 33°682 28-385 2°231 
19°396 36°344 24°839 2-008 
27°114 40°599 17°811 1°45 
357548 43°55° 11-981 0°835 
44°4608 45°516 8-000 o-416 
53714 46°865 5°735 0176 
63°193 47°889 4°651 0°064 
72°860 48°768 4°211 0°020 
82-696 49°592 4059 0°005 
92°695 50°398 47014 ooo! 
102°855 51°200 4°003 "000 
113°175 52°000 : 0000 
123°655 52°800 
134°295 53°600 ; 0*000 
145°095 54°400 , 0°000 
156°055 55°200 ; 0000 
167°175 56-000 “ "000 


NO BAL NO 


Aen 6 & 


6 & 


The form of F for large 7 is 
F = %y*—«"), 
where a’ = 10-602. 


a = 100 R = 10,000 
F F’ Fr’ 


0-000 "000 378-028 
1891 37°816 378°167 
7°562 75°586 376°882 
16°999 113°099 372°786 
30°162 150°008 364°621 
46°966 185°854 351°386 
67°279 220°097 332°491 
gogi2 252°163 307°899 
117°621 281-508 278°204 
147°108 307°675 244°618 
179°039 33°°359 208-837 
248-808 364°974 138°417 
324°157 386°574 80°485 
402-788 398°443 41°431 
483°135 404° 307 19°680 
564°309 407°092 9*605 
645-886 408-556 5°795 
727°700 409°546 4°439 
809°695 410°392 4°096 
891-855 411201 4°017 
974°175 412°003 4°003 
1056-656 412°803 47000 
1139°296 413°603 47000 
1222°097 414°403 4°000 
13057058 415°203 4°000 
1388-178 416°003 4°000 
The form of F for large 7 is 
F = %(n*—a’), 
where a’ = 100-608. 
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SUMMARY 
Further consideration is given to the theory of flow of the most general linear 
elastico-viscous liquid in a coaxial-cylinder elastoviscometer. Useful approximate 
formulae are obtained by employing a method due to Markovitz (1), which is shown 
to be useful in practice only when the annular gap is small. 


1. Introduction 

A NUMBER of instruments have been described in the literature which 
have proved useful in the study of elastico-viscous liquids. Some of them 
(2, 3, 4) are of the coaxial-cylinder type, in which the outer cylinder is 
made to oscillate at a given frequency and the resulting motionof the 
inner cylinder recorded. The theory governing the motion of such instru- 
ments has been considered independently by Oldroyd (5) and Markovitz 
(1). If the cylinders are of length L and radii r, and r, (r, > r,), and if 
the measured quantities are the ratio 6 of the angular amplitudes of the 
two cylinders and the phase-lag c between them, then the equation 
governing the motion of such an instrument is most conveniently ex- 
pressed in the form 


{(22pr? Lj) J_(ar )¥, (arg) — Yq(ar, oJ, (ar,)) + 
ee — 4 {(K 4r?n*) — I} J, (ar, )¥, (a2) — Y,(ar,)J,(ar,))] 


BO se Rind : 1 
6 | (2apr? r, L/a«)(J,(ar,)¥,(ar,)—Y,(ar, J, (or,)) | (1) 





where / is the moment of inertia of the inner cylinder about its axis, K is 
the restoring constant of the torsion wire, p is the density, and a is given by 


x 


fi [ N(r) dr 
—2minp - ——., 
J (1+20nr) 
0 
n being the frequency of oscillation and N(r) the relaxation spectrum (6). 
Equation (1) may be simplified slightly by expressing the Bessel functions 
[Quart. Journ. Mech, and Applied Math., Vol. XIV, Pt. 4, 1961] 
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in the denominator as a Wronskian determinant and using a result found 
in (7). The equation then becomes 


wr: 


te 2 
5 - = [Jp(ar )¥, (arg) —¥a(ar),(are)]+ 


. . 
spac Tagan — 1) Milos) Yorba] (3) 

This equation is unsuitable for computation in its present form since 
there appear to be no tables of J, and Y, for general complex argument. 
Consequently, approximations have to be made to facilitate computation. 
These approximations depend on the size of the Bessel function argu- 
ments, which are accordingly divided into three artificial regions. Large 
arguments are defined as those permitting the use of asymptotic expan- 
sions for the Bessel functions and small arguments as those permitting 
the use of power series for the Bessel functions. Arguments of inter- 
mediate size are then formally defined by Markovitz (1) as being such that 
one can use neither the power series valid for small arguments nor the 
asymptotic expansion which is valid for large arguments. 


2. Approximate formulae 

An approximate formula for large arguments of the Bessel functions 
has been obtained by Oldroyd (5) by substituting for the Bessel functions 
their asymptotic expansions in powers of (ar). In this way, equation (3) 
reduces to 


= a (2)![{1 en we coe x('g—1,)+ 


~) ia 
r 875 





15r,—3r, 15 3 45 os 
a+ ——-— ——,.]S|s a ee 
+| 82 ax * east ("+ aac ed “ses rs ro| (4) 


where ro = (7,72) and S= scram!) (5) 
This equation has been used extensively by Oldroyd, Strawbridge, and 
Toms in the interpretation of their experimental results concerning dilute 
polymer solutions [see, for example, (2, 5)]. 

A similar treatment for small arguments of the Bessel functions is not 
to be found in the literature. However, an approximate formula for this 
range can be obtained without any difficulty by substituting the power 
series expansions for the Bessel functions in equation (3). In this way, 
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the following approximate formula is obtained: 


ee {I—(K/4x®n2)} (r? (r3—r?)? 
gid as (ia)* a |t 


{1—(K/4x*n*)}/, 
+a 32npL («1 





r 
Og. a4 4— 


2 1 





"7 
1+ 
r2 


1 


+ T927 4 


—r)(r§— 5r3 r?— 2rt) + 12r3 rf log, “|, (6 


where terms in a® and all higher terms have been neglected. 

The two cases considered so far present no real difficulties since they 
lend themselves to obvious approximations. The difficult region is clearly 
that between these two extremes, where the arguments of the Bessel 
functions are of intermediate size. Markovitz (1) has attempted to obtain 
an approximate formula for this range by using Taylor series expansions 
and other properties of the Bessel functions. Fundamentally, his analysis 
is equivalent to making a Taylor expansion of the Bessel functions with 
argument ar, about ar,. Equation (3) has first to be recast in the form 


,1¢ 


é my a® K , 4 
ee E + dpr Hass 2) [Elon (or) Yaron) + 


2r, | 4pr,r, L\4n°n 


+a a 


~ Llara)¥ i(07,)—J4(ar,)¥,(arg)], (7) 


where J{(ar,) is the value of a ee at r = r,, and the combinations 
of the Bessel functions occurring in this equation can then be expressed 
as double series in the form 
: T~s—f,\" , 
J,(ar,)¥,(ar,)—J,(ar,)¥,(ar,) = :> (cor, )*4 > ors 3 4 bmajits (8) 
am" 4 
with a similar expression for J,(ar,)Y}(ar,)—J4(ar,)¥,(ar,). Without intro- 
ducing any approximation into the theory, Markovitz is able to evaluate 
these sums and the final equation is expressed as a series in a*. As is to 
be expected, this equation [equation (55) of (1)] is precisely the same as 
equation (6) above. In other words, Markovitz has simply obtained the 
power series formula in a very complicated way. We must conclude, 
therefore, that the intricate mathematical analysis involved in the deriva- 
tion of equation (55) of (1) is of doubtful value. This is not a criticism 
of the fundamental technique of taking Taylor expansions of the Bessel 
functions, but only of Markovitz’s application of this technique in the 
problem under consideration. It will now be shown that the technique 
can be extremely useful when the annular gap is small. 
Equation (6) contains the first three terms of the infinite series (con- 
5092.56 rf 
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vergent for all values of «) which would arise if an exact equation for e/@ 
were required. The rapidity of convergence of this series will therefore 
determine the usefulness of equation (6) as a reliable approximation. It 
is clear that the series will converge rapidly for small «, so that equation 
(6) is a reliable approximation for small arguments of the Bessel functions. 
It also appears likely that the series will converge rapidly for fairly large 
values of a provided that (r,—r,) is small enough. In other words, it is 
possible that equation (6) will be a reliable approximation for any value 
of the Bessel function argument provided that a(r,—r,) is small. This 
suggests that we follow Markovitz’s method of expanding the Bessel 
functions as Taylor series, taking «(r,—r,) as the small quantity and not 
xr, and ar, as in the derivation of equation (6). In practice, this means 
that in the series 


J,(a'g) = > lata Pa pear.) (9) 
m=0 ~ 


and Y,(ar,) = > fet Y{"(ar,), (10) 
m=0 y 


only the first few terms are to be considered. The next section illustrates 
the ease by which approximate formulae can be obtained when the annular 
gap is small enough to allow the Taylor series to be terminated in this 
way. 


3. Special case 


As an example, we consider the case when the annular ,ap is small 
enough to allow the series to be terminated at m = 2. In this case, the 
combinations of the Bessel functions occurring in equatio:. (7) can be 
written in the form 


J,(ar, )¥, (arg) —Y,(ar,)J,(arg) = a(r.—r,)[J,(ar,)V4(ar,)—Y,(ar,)J4(ar,)|+ 
+ $07(r.—r,)"[ (ar) ¥i(or,)—Ky(ary)Jy(ary)}, (11) 
J,(arg)¥ (ary) —J4(or Yj (arg) = [J,(ary)¥y(ar,)—J4 (ary )¥y(ar,)]+ 
+-$a2(r.—r,)*{ Jj(ar,)¥4(ar,)—J4(ar,)¥y(ar,)|, (12) 


and using results found in (7), these reduce immediately to 


J,(ary)¥, (ar) —Y,(ar,)J,(ar,) = Feene _ ). (13) 


2 
mr; mrt 


‘ 9,8 
and = J, (ar,)¥j(ar,)—J}j(ar,)¥,(ar,) = |e al (14) 


wor; mr, ar 
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Equations (13) and (14) can now be substituted into equation (7), and 
after some rearrangement, the following approximate formula is obtained 





. a pe] A Pt) ta. (15) 


3 5) 
4npr, 1, L ry ri 2r, 


This equation could be used to interpret the experimental results 
obtained from the elastoviscometer of Oldroyd, Strawbridge, and Toms 
(2), since one of the annular gaps available in this apparatus is as small 
as 0-05 cm and for the liquids investigated by Oldroyd, Strawbridge, and 
Toms, |a| < 4. Even with an annular gap of 0-1 em, which is also available 
in the apparatus, equation (15) could be used as a reasonable approxima- 
tion over the major (and most important) part of the frequency range. 

The use of equation (15) would greatly facilitate the interpretation of 
the experimental results of Oldroyd, Strawbridge, and Toms, since these 
have been interpreted hitherto by means of the complicated asymptotic 
expansion formula [equation (4)]. However, since amplitude ratio alone 
was measured in their experiments, it would be necessary to take moduli 
of both sides of equation (15) before the equation could be used to interpret 
the experimental results, and, in consequence, it would not be possible to 
derive « as a direct function of the frequency n. This means that the 
interpretation of the experimental results would again have to follow the 
laborious procedure of the preceding papers (5, 6) of comparing families 
of theoretical (@,”) curves for certain idealized spectra with the experi- 
mental curve. The computation involved in this comparison would ob- 
viously be simplified by employing equation (15) in place of equation (4), 
but this is the only improvement which can be expected when amplitude 
ratio alone is measured. 

On the other hand, if both amplitude ratio and phase-lag results were 
available, the use of equation (15) would introduce a significant improve- 
ment into the interpretation of the experimental results. In this case, it 
would be possible to derive a as a direct function of the frequency n, and 
the relaxation spectrum N(r) could then be obtained from equation (2). 
This adds further weight to a conclusion reached in Paper II of this series 
(6) that measurements of phase-lag are indispensable if a complete charac- 
terization of elastico-viscous liquids is required. 


4. Conclusion 

As a general conclusion, we may say that the technique of expanding 
the Bessel functions as Taylor series is of no practical importance when 
all the terms in the series are considered, since the resulting formula can 
be obtained in a much simpler way by substituting the power series 
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expansions for the Bessel functions in equation (3). It will be useful in 
practice only when the annular gap is small (i.e. when the Taylor series 
can be cut short at some small value of m). In this case, it eliminates 
the necessity of considering the rapidity of convergence of the infinite 
series of which equation (6) contains the first three terms. The value of 
m chosen will depend on the apparatus under consideration, and quite 
obviously the complexity of the resulting equation will depend on the 
value chosen for m. For example, when the series is cut short at m = 3. 
the approximate formula obtained is of the form 

Or 1. fe Ree (r2—7,)*\ __ rte a) a 
6 4npr,r, L ry r r } 2r, 





af{(K/42?n?)— 7} 

l2mprir, L 
which is clearly more complicated than equation (15), and the interpreta- 
tion of experimental results would be made correspondingly more difficult. 





(r2—1,)°, (16) 
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SUMMARY 

[t is supposed that a localized distribution of current is created in a viscous, 
conducting, incompressible fluid initially at rest in a uniform magnetic field. 
' A mechanical disturbance is applied simultaneously, and the resulting motion is 
then allowed to decay. The asymptotic forms of the velocity and magnetic fields 
are obtained, and it is found that in the final stages of the decay the velocity field 
is in general due to two ‘viscous vortex rings’ moving with the Alfvén wave velocity 
in opposite directions along the undisturbed lines of force. The rings are not 
identical, their axes being inclined to one another and their intensities unequal. 
The difference in strength depends only on the initial distribution of current. With 
each ring is associated a current distribution decaying at the same rate as the 
velocity field. The current lines for one ring are parallel, and for the other anti- 
parallel to the vortex lines of the ring. 


1. Introduction 

Ix a recent paper, Saffman (1) has extended the work of Phillips (2) on 
the decay of a general localized disturbance applied to a viscous incom- 
pressible fluid initially at rest to the case in which the fluid is an 
electrical conductor permeated by a uniform magnetic field. In this 
extension, it was supposed again that the disturbance was mechanical, so 
that its effect was initially to impose a velocity field of limited extent on 
the fluid without disturbing the magnetic field. That is, the disturbance 
was regarded as giving rise to a localized vorticity distribution without 
producing any current at the initial instant. Subsequently current is in- 
duced because of the interaction of the uniform magnetic field with the 
velocity field. 

In the present note we extend this work further to the case when the 
disturbance has electromagnetic effects as well as mechanical ones, so 
that a localized distribution of current, or current blob, is created initially. 
As previously, we shall examine the form of the velocity field in the final 
stages of decay, and it is now found that the results are modified. Two 
decaying ‘viscous vortex rings’ are again produced, which move in opposite 
directions along the undisturbed lines of force with the Alfvén wave 
velocity, but they are no longer symmetrical. Their axes are now not 
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parallel and each ring has different amounts of linear momentum, the 
difference being proportional to the volume integral throughout space of 
the magnetic field perturbation. 


2. B-haviour of the disturbance at infinity 

W:. employ the same units and notation as in (1) to which the reader 
is rearred for details. It is supposed that at the initial instant the per- 
turbation in the magnetic field is due to a localized distribution of current, 
which must be such that 


i) jdx =0 
since the current lines will be closed. It now follows from the expression 
for the magnetic vector potential in terms of the current distribution that 
initially 
- 
where H is the magnetic field and H, its undisturbed value. 


Let us now consider the order of magnitude of the current field at in- 
finity. The equation for the current is 


h = O(r-*) asr>o, (1) 


a = S curl curl(u A H)+AV4¥j, 
ot 20 4n 
where A is the magnetic diffusivity. The second term on the right-hand 
side of (2) represents the diffusion of current due to the conductivity of 
the fluid, which by itself produces an exponentially decaying distribution. 
The first term on the right-hand side describes the production of current 
due to the interaction of the velocity and magnetic fields. Since [see (1)] 
u = O(r-*) and H ~ H, as r > ~, it follows that aj/ét, and hence j, are 
at most O(r-) at infinity. This current distribution at infinity will produce 
a magnetic field of order r~* at infinity, so that (1) remains true at subse- 
quent times. 

The argument given in (1) for the order of magnitude of the vorticity w 
at infinity still applies, so that w = O(r-*)t. It follows that the expansion 
obtained previously for the Fourier transform of the velocity field, namely, 


p(x, t) = (2m)-* [ u(x, the-™-* dx, 


is still valid, and we have 


K 


b,(x,t) = (5) —“) 5+ 0tW, (3) 


where M = (27)-* { udx. We denote the Fourier transform of the current 


+ In fact, a more careful analysis shows that j = O(r-*) and w = O(r-*), but the orders 
of magnitude given above are sufficient for our purposes. 
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by (x, ¢). Then in virtue of the behaviour of j at infinity, the first deriva- 
tive of ~ at « = 0 exists and, since divj = 0, we have the following ex- 
pansion [see (2)}: iH 
$; = Fy Sue's £,.+0(x), 


where € may depend upon the time ¢ but is independent of x. 

From the relation between h and j, it follows, by arguments identical 
with those in (2), that I(x,¢), the Fourier transform of h(x,?t), has the 
expansion 


i= (5 —“P) +010 (4) 


K 


about « = 0, and, moreover, that 


PS» a 
t= | h dx. (5) 


The non-analytic behaviour of T at« = 0 reflects the fact that the integral 
(5) is conditionally, but not absolutely, convergent. 


3. The final stage of decay 
The equations of motion can be reduced to a pair of non-linear simul- 
taneous equations involving m and I [equations (13) and (14) of (1)]. 
Expanding the terms in these equations about « = 0, we find, as in (1), 
that the zero-order terms give 
i he on 
dt dt 


That is, M and €¢ are constants of the motion whose values are determined 
by the initial conditions. The physical significance of M is immediate, 
since 87%pM is the net linear momentum of the fluid (p is the density). 
That of ¢ is not yet clear, although, as will be seen, it is associated with 
the spatial distribution of momentum in the final stage of decay. 

After a sufficiently large time, the amplitude of the disturbance becomes 
smal! because of the dissipative action of viscosity and conductivity, and 
we may then neglect the non-linear terms in the equations of motion. 
These, or rather their Fourier transforms, then reduce to 

Fe. nth, —iat,T, = 0, Tb actly — ing gy = 0, 

ot ot 
where the l-axis has been taken in the direction of the undisturbed magnetic 
field, a? = »H?/4xp is the square of the Alfvén velocity, and v is the 
kinematic viscosity. These equations can now be solved readily to give 
ep and IT at subsequent times in terms of their values at some instant t = ft, 
for which the non-linea’ terms have become negligible, as in (1). 
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The velocity and magnetic fields are given now by the inverse transform 
relations as integrals over wave number space involving @(x,¢,) and 
I'(x,¢,). Now the integrands also contain a factor exp{— }«*(A+-v)(t—t,)}, 
so that when ¢—t, is large only the values of e(x,t,) and I(x,¢,) near 
« = 0 contribute to the integrals. Substituting the leading terms of the 
expansions (3) and (4) for ~(x,?,) and I(x, ¢,), we find [as in (1)] that the 
asymptotic forms of the velocity and magnetic fields as ¢t + 00 are given 
by the Fourier transforms of the following e eaggas 


(x,t) = (5—“ saa | Y <2 sinh 107] M,+ 


—— = 


—— (; sinh 407, (6) 


K 


T(x, t) = (8,,—““d]e-4*0+| [ cosh apa Or)l,+ 
ar. CS) : 


2, 
+288) M, sinh 107], (7) 


where + = t—t, and © = {x*(A—v)?—4a%«j}. 

The asymptotic forms of the velocity and magnetic fields in the final 
stage of decay thus depend upon the initial conditions only in so far as 
the latter determine M and ¢, and are independent of the history of the 
decay. 


4. The velocity and magnetic fields 
When the two diffusivities are equal, i.e. A = v, the expressions (6) and 
(7) simplify considerably and the Fourier transforms can be evaluated in 
closed form. We then find, exactly as in (1), that the vorticity distribution 
in the final stage of decay is 
w(x, t) = —(M-+af) a V®,—(M—aBd) A V9,, (8) 


where 


3 
®, = (=)Fexp{—lt@ , tar)? +23+-23]/47 v4) 


‘ (9) 

3 

®, - 1(=)Pexpt —[(x,—ar)*+ 23+ -23]/4v7} | 
VT 


Similarly, the current distribution is sae by 
j(x,t) = — oo M+ at) A Vo, +2 — (M— at) A V,. (10) 


The contribution to the velocity field from the term involving ®, is 
that of a viscous vortex ring, of the type described by Phillips, whose 
centre is moving anti-parallel to the undisturbed magnetic field with the 
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Alfvén velocity and which decays as if it were in a non-conducting fluid 
of kinematic viscosity v. The linear momentum associated with this ring 
is }(87°p)(M+<af). The axis of the ring is parallel to its limear momentum, 
and the vortex lines are circles centred on the axis. The velocity field 
arising from the term involving ®, is likewise a viscous vortex ring, but 
its centre moves in the opposite direction and the associated linear 
momentum is $(87%p)(M—af). 

When the initial disturbance is purely mechanical, as in (1), § = 0 and 
the two rings are identical, provided M + 0. If the disturbance is both 
mechanical and electromagnetic, then the two rings are not identical, 
their axes are not parallel, and their intensities are unequal. The initial 
current distribution is reflected in unequal amounts of linear momentum 
for the two rings. When M == 0, which is the case if the disturbance is 
purely electromagnetic or the net linear momentum is zero, then the rings 
are mirror images of one another. If the disturbance is such that both M 
and © are zero, then it can be shown that the asymptotic form depends 
upon the entire history of the motion and a universal asymptotic state 
does not exist. 


Comparison of (8) and (10) shows that the current lines associated with 
each ring are respectively parallel and antiparallel to the vortex lines, aid 
similarly for the streamlines and the lines of force of the magnetic field 


perturbation. 

For the general case when the two diffusivities are unequal, the analysis 
is more involved, and the Fourier transforms of (6) and (7) cannot be found 
in closed form. However, approximate methods of integration can be 
employed and it can be shown (the analysis is straightforward but cum- 
bersome) that the same conclusions still hold provided a?r/ |A—v| > 1, 
except that now the rings decay as if the diffusivities were both equal to 
their arithmetic mean. This last condition will always be satisfied even- 
tually, and then the expressions (8), (9), and (10) will be valid with v 
replaced by $(A+-y). 
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SUMMARY 


It is shown that a class of axially symmetric boundary value problems in potential 
theory for the region external to a lens may be reduced to problems in two-dimen- 
sional potential theory. The solution for the general axially symmetric Dirichlet 
problem is obtained and it is shown that the Stokes stream function for an arbitrary 
flow may also be reduced to a Dirichlet problem in potential theory. The Green's 
function for the lens is also derived and the extension of the approach to non-axially 
symmetric problems is considered. 


1. Introduction 

Ir has recently been shown (1, 2, 3) that a class of problems for Helm- 
holtz’s equation with boundary conditions applied on a wedge can, by 
employing a suitable integral representation for the solution, be reduced 
to the solution of a problem in two-dimensional potential theory. In view 
of the powerful tools available in two-dimensional potential theory it 
seems of interest to examine whether there exist other classes of boundary 
value problems which can also be so reduced. In the present paper a class 
of axially symmetric boundary value problems for Laplace’s equation is 
considered from this point of view. The problems treated are those where 
the boundary conditions are to be applied on a lens. 

It is shown that, by employing a suitable integral representation for 
the solution, the general axially symmetric Dirichlet problem for Laplace’s 
equation in the region external to the lens is reduced to a simple two- 
dimensional Dirichlet problem. The general approach is developed in 
section 2 and is applied to determine the potential function which is 
constant everywhere on the lens surface. The method is also extended 
to solve the general Dirichlet problem for the lens. It is shown in section 3 
that the solution of this latter problem may be employed to obtain the 
Stokes stream function for arbitrary axially symmetric flow past the lens. 
The direct relationship obtained between the stream function and a poten- 
tial function seems to be new though a relationship between the stream 
function and a potential function in five-dimensional space has been used 
with great success by Weinstein (4) and Payne (5). 
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In section 4 the general approach of section 2 is applied to determine 
the velocity potential for a uniform stream flowing parallel to the axis 
of symmetry past the lens. The application of the method is straight- 
forward but the manipulations necessary are more complicated than those 
for the Dirichlet problem. The final solution obtained agrees with the one 
obtained by Shiffman and Spencer (6) who applied a difficult and in- 
genious procedure involving the method of images in a multiply-sheeted 
Riemann—Sommerfeld space. 

The general method of reducing the three-dimensional potential problem 
to a two-dimensional one is not confined to axialiy symmetric problems, 
and in section 5 it is applied to determine the Green’s function for the 
lens. In principle this Green’s function could be used to solve the Dirichlet 
problem of section 2, but the direct approach seems to be simpler. It is 
shown in section 6 that the solution of the non-axially symmetric Dirichlet 
problem for the lens can be expressed in terms of a sum of solutions of 
the axially symmetric Dirichlet problem. This result is applied to the 
particular problem of a conducting lens placed in a uniform electrostatic 
field which is perpendicular to the axis of symmetry. It is shown that the 
solution may be expressed in terms of the solution to the first boundary 
value problem solved in section 2. For the particular case of a circular disk 
the solution reduces to one given by Hobson (7). 

Problems of the type treated in sections 2, 3 have also been considered 
by Collins (8); he also solves the problem by employing a convenient 
integral representation for the solution, but the approach is different from 
that presented here. 


2. Dirichlet problems for a lens 

The most convenient method of specifying a lens is in terms of toroidal 
(or peripolar) coordinates o, %, d. These are related to cylindrical polar 
coordinates p, ¢, z by the equations 


asinho asin ys 
~~ cosh e—cos x’ 


~ 
~ 


p= 


cosh o—cos yb’ 


where a is constant. The surfaces 4 = constant are spherical caps and 
we shall consider the lens formed by the intersection of the caps y = a—f 
and % = —B, where a > 8, 8 > 0. For 8 < w both surfaces of the lens 
will be convex but for 8 > 7 one surface will be concave. The region 
external to the lens is defined by —B <% <a—f,0 <a <ao;0=0 
(y% ~ 0) corresponds to points on the z-axis and o = 0, % = 0 refers to 
points at infinity. At points on the rim of the lens o becomes infinite. 
We shall first consider the determination of the electrostatic potential V 
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outside the lens when the surfaces are maintained at a constant poten- 
tial U. In this case V is independent of ¢ and satisfies the equation 


sinh o é} é (sinh o eV 
= ae a =) “— () 


where s = cosho, tr = cosy. V must be finite everywhere at points out- 
side the lens. It must also tend to zero at infinity and be equal to U on 
% = —B and ¢ = a—f. 
Writing V in the form eae then 
ew 


— W = 0. 
+ +55 


It can be verified directly that a solution of (2) is 


— ( g(w, yp) de 


~~ J (cosh w—cosh a)!’ 
oc 





(3) 


e*g eg 
ew? * af? 
It is then easily shown that equation (3) represents a possible solution of 
(2) provided that g is bounded as w — 00 and is non-singular in the region 
—B <% <a—B, w >0. In order that (s—r)!W, where W is defined by 
equation (3), should be non-singular on the axis of symmetry (o = 0) we 
require that the integral for W be finite when o = 0. This implies that 
g(0,%) = 0 which also ensures that V +0 at infinity. The boundary 
conditions on the lens give 


[ gw,—A)dw 
J [coshw—cosho]! [cosh e—cos f}*’ 


where = 0. (4) 





g(w,a—B)dw RS: 
[cosh w—cosho]}* [cosh o—cos(a—8) |” 





Equations (5) and (6) are of the Abel type and may be solved to give 
—B) = 


sinhw (7) 
= (aaa w—cos f}’ 
By! sinh w 
a [cosh w—cos(a—)] 
All the conditions of the original boundary value problem will thus be 
satisfied by determining a particular solution g(w, /) of equation (4), finite 
everywhere in 0 << w <0, —B <4 <a—fP. g(w,y) also satisfies equa- 
tions (7) and (8) and vanishes when w = 0. This boundary value problem 





g(w,a—B) = (8) 
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for g is of a well-known type. Its solution is easily obtained by classical 
methods and is given by 


@ 


gw.9) = Meh f|____ fio 


cosh n(w—t)—cos n(y-+ 8) 


f(t) | dt, (9) 


cosh n(w—t)+ cos n(y+-B)} 





—-@ 





where n = m/« and f, and f, denote the right-hand sides of equations 
(7) and (8) respectively. The integral for g(w,) may be evaluated to give 


£f[nw, n(yb-+28)]—f[neo, mi} +o f{w, #], (10) 


vy 
g= — 
x 
sinh w 
r} = Dp. =  ——_.. 
eer fiw. #] cosh w— cos x 
If we substitute g from equation (10) into equation (3), it can be shown 
that 





ie f pe tk * {Fi nw, n(yb-+28)]—f[nw, nf} dw) 
.=— .* ree | (cosh w—cosh o)! \ an) 


oa 


This form for V agrees with a result obtained by Macdonald (9). 

The above analysis may be extended to obtain the solution of (1) when 
the values of V prescribed on the lens are arbitrary functions of o but are 
independent of the angular variable ¢. The other conditions imposed on 
V are still assumed to hold. If V(o,—8) = U,(c) and V(o,a—f) = U,(c) 
are the prescribed values, then equations (5) and (6) will again hold if 
in the right-hand sides U is replaced by U,(e) and U,(c) respectively. 
The right-hand sides of equations (7) and (8) will now be replaced by /,(w) 
and f,(w) respectively, where 


2D 
” 


U,(c)tanh o do ) 
| cosh o—cos 8 }*| cosh o—cosh w}}} ‘ 


as! (cosh w)* d { 
f(wv) = -— ae aaa w)t 





(12) 
U,(c)tanh o do ). 

[ cosh o—cos(a—)|*{ cosh o—cosh w}}| 
(13) 


__ (cosh w)! d { ; 2 
f,(w) = a ae ibe uw) | 





Since U, and U, are only given for o > 0, equations (12) and (13) define f, 
and f, only for w > 0. It can be shown, however, that /,(0) = f,(0) = 0, 
and f,(w) may thus be defined for w < 0 as —f,(—w); a similar definition 
holds for f,(w). The solution g will once more be given by equation (9). 
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3. The stream function for a lens 
The equation satisfied by Stokes’s stream function S in cylindrical polar 


coordinates is 2S 1eS as 
Sep een aetboemen ann i (14) 
ep é2* 


The general mathematical problem of determining S for a lens is that of 
obtaining a solution of equation (4) bounded everywhere at points ex- 
terior to the lens which tends to zero at infinity and takes prescribed values 
on the lens surfaces. We shall assume that S = S,(c) on ¢ = —f and 
S = S,(c) on & = a—B, where S,(0) = S8,(0) = 0. This assumption leads to 
no loss of generality (8). It may be verified that S may be written as 
(cosh o—cos ys)? sinh o(@W /éc) where W satisfies equation (2). Thus it is 
seen that all the conditions of the problem may be satisfied if a solution 
W of equation (2) can be obtained satisfying all the conditions imposed 
on W in section 2 and such that 
W(o, —8) = [ S,(w)| cosh w—cos B}'cosech w dw, (15) 
W(e,a—B) = [ S,(w)| cosh w—cos(a—)]'cosech w dw. (16) 


a 


The boundary value problem for W is of the same type as that for the 
potential problem of the previous section. Thus the stream function may 
be expressed in terms of the solution of a Dirichlet problem for equation(1). 

For the particular case of a stream of uniform velocity U flowing in the 
z-direction past the lens, S, and S, are given by 


, 1  Ua'sinh*o “ae Ua? sinh?o 
ac ae [cosh o—cos B |?’ - § [cosh o—cos(a—B) |?” 





The right-hand sides of equations (15) and (16) now become —a?x the 
right-hand sides of equations (5) and (6), and thus 


S = —a*(s—r)- sinh o—— {V(s—r) }. (17) 
Co 


where V is given by equation (11). This agrees with the result obtained 
by Payne (5). 


4. The ve‘ocity potential for uniform flow past a lens 

In this section the velocity potential V is obtained for the case of a 
uniform stream of velocity U at infinity and in the z-direction flowing 
past the lens. V must satisfy all the conditions imposed on the potential 
of section 2 except those on the lens, where now @V /@b = U éz/ép. 
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The integral representation of section 2 will again be employed so that 


5 = (s—r)* | sa (cosh w—cosh 0) dw +- 


+$siny(s—r)- { g(cosh w—coshc)-+ dw. (18) 


It is convenient at this stage to note that the general solution g(w,) of 
equation (4) may be written as F(w-+i¢)+G(w—i) where F and G are 
arbitrary functions. Partial integration of the first integral in equation 
(18) leads after some simplification to 


7 (cosh w—cos x) 
(cosh w—cosh o)* 





| cosech (w+ iy = {sinh }(w-+-ib) F(w+ub)} — 


—cosech }(w— ip) “ {sinh }(w—tp)G(w—wb)}— 


sinh w{ F(co)— ee) 
ban , (19 
faa i 
where F(oo) and G(co) denote the values of F and G when w+. The 
analysis leading to equation (19) may be shown to be valid provided that 
F and G and their first derivatives are bounded and continuous and that 
F and G satisfy a condition of the form 


|K(w)—K(o0)| _ 
jwo|—> 20 exp{—(e+4)|\w}} 


where C is bounded, K stands for either F or G, and « > 0. 
If F, and G, are defined by 








(20) 


F(w) = sinh joA (F, cosech $w). G(w) = sinh joo (G, cosech 42), 
then it is easily verified that equation (19) becomes 


ov ; ° (cosh w—cos ¥) 
——s one —-7)t 
Oy a iii (cosh w—cosh a)? 





| DER C0-+ig) — G,(w—ig)}+ 


sinh w{F,(00) — G,(00)} 
4(cosh w—cos ys) 





| dw (21) 


where D is the operator d?/dw*—}. 
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The boundary conditions on the lens give two integral equations which 
may be solved tégive 
D{F,{w—i|—G,[w+ip} 
_ Ua # sinh w 1 sinh w{ F,(00) —G,(%)] 
~~ qi dw*[coshw—cosB] 4  [coshw—cosf] 
Di F[w+i(a—Bp)|—G,[w—i(a—B)} 

_Uad __sinhw at sinh w| F,(00) —G,(o) | 
~~ qt dw* {cosh w—cos(a—f)| 4 [cosh w—cos(a—B)] © 
Solutions of equations (22) and (23) are 
F[w—ip]—G,{w-+ ip} 

Uaj_ sinhw 


a — k(w, —p)| +4[F,(00)—G@,(20) Jk(w, —B), (24) 


mi \cosh w—cosB 


F.[w+i(a—)]—G,w—i(a—B)] 
hed ‘aj sinh w woe iA 
= mi \cosh w —cos(a— wen B)\ + 


+4[F,(00)—G,(00)]k(w,a—B), (25) 





» (22) 





(23) 


where 


k(w, ys) = sinh 4w cos $f log( oe $+ cos ) a 





cosh 4w—cos $y 
+sin }% cosh $wtan~(sin 44 cosech }w). (26) 
The right-hand sides of equations (25) and (26) are particular integrals 
of equations (22) and (23) subject to the condition (20). The integrand for V 
must vanish when w = 0; a sufficient condition for this is F,(w) = G,(—w). 
The latter condition will be satisfied if F,(0o0)—G,(0%o) = Ua/mi, and in 
this case 


F[w—ip|—G,|[w+ip] = {gah aoa B)’ 


Ua sinh w 
mi [cosh w—cos(a—B}]" 


F[w+i(a—B)|—G,[w—i(a—Bp)] = (28) 





The problem thus reduces to determining a potential function 


g(w, f) = F(w+ih)—G,(w—ip) 
vanishing on w = 0 and with values on ¢ = —f and % = a—f given by 
equations (27) and (28) respectively. This problem was solved in section 2 
and hence 
F,(w+i)—G,(w—ig) = = {flme, (y+ 28) |—f[nw, mp} hae 22 pw, $), 
(29) 
Gg 
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where f[w,#] is defined in section 2. Hence 


Fw) = = coth yw +o {eoth 4n(w+ 21B)—coth jnw}, (30) 


G,(w) = F,(—w). (31) 
If the values of F,(w) and G@,(w) given by equations (30) and (31) are 
substituted into the equatjon for V we obtain 
V = Uz+(s—r)! [ [sinh g(wo-+ ip) {UP igeosech (+ if)} + 


~ 


‘ ; dw 
+-sinh }(w — i) <(F (ub—w)cosech }(w— iw com waa 998 
(32) 
where F*(w) is the second term on the right-hand side of equation (30). 

Equation (32) can, after some simplification, be shown to agree with 
the contour integral form of the solution given by Shiffman and Spencer 
(6). 

In the above problem the form of éV/@ is such that F, and G, can be 
found fairly easily. For arbitrary values of @V/é the necessary analysis 
is not readily obtained. The above approach may be used to reduce the 
general Neumann problem to the solution of equations of the form (22) 
and (23). It does not seem possible to obtain a general method of deter- 
mining the particular integral of these equations such that the integrand 
for V is zero at w = 0 and F, and G, satisfy equation (20). It appears that 
this part of the solution has to be carried out separately for each problem. 


5. The Green’s function for a lens 

In this section we shall derive the Green’s function G for the lens. The 
general boundary value problem of section 2 could be solved in terms of 
G but in view of the involved integrations in Green’s formula the Green’s 
function method is far less direct than that of section 2. The determina- 
tion of the Green’s function is presented here as an application of the 
present approach to a non-axially symmetric problem. G is defined as 
G,+1/R, where R represents the distance from a fixed point Q and G@, is 
a solution of Laplace’s equation, i.e. 

eam a é oy 52 fe “4 4, 1 eG, a (33) 
éo\ s—r 0) ' @b\s—r ab} * (s—*)sinho ad? 

G, must be bounded everywhere, and tend to zero at infinity and @ must 
vanish on the lens. Physically, @ represents the electrostatic potential 
due to a point charge at Q when the lens is earthed. 
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If oo, uo, dy denote the coordinates of Q then 


1 1 (cosha—cos #)*(cosh o)—cos yg) 
R- av2 [cosh 8—cos(xs—yp) |* 
where cosh § = cosh o cosh o,—sinh o sinh o, cos(¢— dy). 


Also, if G, = bier. then 





(34) 


ee? PG, , &G, 
Gs tho G. tt he? ae @. 35 
Since ¢ me occurs in the cosh 6 term in annie (34) it seems reasonable 
to consider a solution of (35) depending only on ¥% and 8. It can then be 
shown that G,(5,%) satisfies 
2G. é "Gi, _ 
57 +00 coth 3 440,45 aad (36) 


A solution of equation (36) is 


Q , g(w, f) dw 
13> = 
(cosh w—cosh )*’ 


where g is a solution of equation (4). In order that G, should be bounded 
everywhere we require that 9(0,) = 0. 
The boundary conditions on the lens give 
(cosh o,—cos yy)! sinh w 
(w, —B) on r i’ 
ann 2| cosh w—cos(B +-y,) | 
er “pie (cosh o,— cos #,)' sinh w 
~ — anv2[ cosh w—cos(a—B—yo) | 





(37) 





(38) 


The boundary value problem for g is precisely that of section 2 and g(w, #) 
will be given by equation (9) if f, and f, are taken to denote the right-hand 
sides of equations (37) and (38) respectively. The integrals may be 
evaluated to give 


g(w, yp) = —(C8h.co— 008 Ho)" ¢ pra, s1(yh-+- o-+28)]—f [n0e, (yb —e) }— 


aan? 


__ (cosh o,—co8 yo)! * fw 


7an2— 


oH). (39) 
Finally 
(cosh o,—cos yy)#(s—r)* | 

aav2 


x f (de . Mb+-o+ 28) |—f[nw, n(y— wp) | py 


| cosh w— cosh 5 }* 


G= — 








(40) 
) 

For the particular case of the spherical bowl equation (40) may be identi- 

fied with a solution obtained by Hobson (7). 
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6. Non-axially symmetric problems 

The analysis of section 2 may be extended to cover the case when on 
the lens V is specified as an arbitrary function of both o and ¢. If the 
values of V prescribed on the lens are to be single-valued functions of 
position then they can be expressed as a sum of terms of the form 
F(a) cosn(¢+¢,), where n is a positive integer or zero. Thus, from 
equation (35), the general problem may be solved if we can solve the 
Dirichlet problem for 


faq tooth oo + U—n* cosech* aU + = 0, (41) 

It may be verified that a solution of equation (41) is sinh"o(é" W/és"), 
where W is a solution of equation (2). Thus, by repeated integration, the 
solution of the problem may be reduced to a boundary value problem for 
W identical with that solved in section 2. The limits of integration in 
the repeated integrals must be chosen so that the solution U is bounded 
at the rim of the lens. The condition at the rim will normally be satisfied 
by choosing s = © as the lower limit in each integral. 

The general approach will be illustrated by considering the particular 


case when y _ Yasinhocos ¢(cosh o—cos B)-* ony = —B 


and V = Uasinhocos ¢{cosh e—cos(a—B)|-!_ on # = a—B. 
Physically this corresponds to determining the potential of an earthed 
lens placed in a uniform field perpendicular to the axis of symmetry. 
The boundary conditions show that the boundary value problem for 
W is precisely that of section 2 if U is replaced by —2Ua. Thus the 
potential is 


V, = —2a cos $(s—r)?2 {V(e—r)-4}, (42) 
ca 


where V is defined by equation (11). The author has been unable to 
discover the solution to the above problem in the literature but for the 
particular case of a circular disk the potential as given by equation (42) 
agrees with that given by Hobson (7). 
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SUMMARY 


This paper is concerned with the stress distribution at a rigid circular disk on an 
infinite elastic layer over a rigid toundation when the disk is subjected to either 
a torque about its axis or a moment about a diameter. In the analysis, which is 
within the classical theory of elasticity, dual integral equations are formed and 
solved using the method of J. C. Cooke (1). Some stress curves are plotted for the 
torsion case. 


1. Introduction 


Tue two static problems to be treated stem from the aualogous dynamic 
problems consisting of the torsional and bending modes of forced vibra- 
tions of a rigid disk on an infinite elastic stratum (2,3). In order to solve 
the dynamic problems the contact stresses, as an approximation, were 
assumed to be those of the corresponding static cases for a semi-infinite 
elastic medium with a periodic time factor. The following investigation 
of the static cases is primarily concerned with finding the actual contact 
stresses and the depths of strata at which these stresses may be assumed 
to differ insignificantly from those corresponding to a semi-infinite elastic 
medium. In each case dual integral equations are formed and solved 
using the method of J. C. Cooke (1). 


2. Torque applied to the disk 


It can be shown (4) that only the tangential component of displacement 
ug(r, 2) is different from zero and satisfies the partial differential equation 
ajl @ ug 
a ru — == 0. l 
ar\ ror aj+ é2? ”) 


Consequently, the only components of stress are 


Oflu ou 
To = pr- ( ) and 1,9 =p—. 
er\ r 


~ 
~ 


The cylindrical coordinates are shown in Figs. | and 2 and yp is the 
modulus of rigidity of the material. 
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To solve equation (1) introduce the Hankel transform 


U(é,z) = | ugrd(€r) dr 


so that 


w= f ve neneras, — ro—u {PED erenae. 


0 0 
From (1) and (2) it is readily shown (5a) that 
U = A(£je*+ Biéjef. 
The boundary conditions of the problem are 
u(r,0)=ar (O<r<a), rtyal(r,0)=0 (a<r),  up(r,h) = 0, 
where a is the radius of the disk, A the depth of the elastic layer, ea . 
is the small rotation of the disk. 


LZZZ2 


oO 





elastic 





h 
MLM hhhhhh MAA 
rigi 





z 
Fie. 1 


To satisfy the boundary conditions (3), B(é) = —A(é)e~** and A(E) is 
determined by the dual integral equations 


[ A(é)(1—e-***)E, (Er) dé = ar (0 <r <a), 
0 
[ A(\(L+e-*™) EU, (Er)dg =O (a <r). (4) 
0 
On setting u = af, A(E)(1-+e-*")€ = aag(u), r = pa and h = fa, equations 
(4) become 


x 


[ g(u) u- tanh BuJ,(pu)du = p (0 <p <1), 


0 
o 


( g(u)J,(pu) de =0 (p> 1). 


0 





° 
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Equations (5) are in a form suitable for a direct application of J. C. 
Cooke’s method (1) for finding g(u). Therefore, let 


1 
g(u) = *u { seesin wea (6) 


With representation (6) the second of the integral equations (5) is satisfied 
and substitution into the first leads to the Fredholm integral equation 


1 
fia), | K(p)f(t)dt = «x. (7) 
—1 


1 
The kernel in (7) is K(p) = | (1—tanhv)cos pu dv, where v = Bu and 


0 
p = (x—t)/8. In deriving (7), use has been made of the fact that f(z) can 
be shown to be an odd function of z. 
The method chosen to find f(x) is to replace (7) by a system of linear 
algebraic equations (6) using equidistant data. The function f(x) is then 
obtained in the form 


M 
f(x) 2. 2 Famer em (8) 


the numerical coefficients a,,,,, depending on the value of 8. Interpola- 
tion (8) by powers is performed by means of Chebyshev polynomials (7). 
Representations (8) and (6) give 


4 mM 
g(u) = ye 2, Dom +1 Lom+1(™); 


u “ OJ, »4g(U). 


where Tem 44(4) = (47m! = rr ae 


To find the contact stress 7,4(r, 0) use is made of the result (5b) 
x p(l—p?)n 
| wi" -DS) .a(u) J (pu) du = 2"-1T(n+4) 


. 0 (p > 1). 
This gives 


T@(7,0) = —pao ( g(u)J,(pu) du 
0 


M M 
4 2"n! m! 
ec ted — 1)"(1 —p2)"-? ——_ 
she 2 y"(1—p*) (2n)! < ~ (m— —nyi mth 
where p = r/a < 1. The torque required to maintain the rotation a is 
M M 
2"(n+-1)! m! 
T = 32yaa* —1)* al sO : 
-~ 2. ) (2n+-3)!  (m—nyim 


Numerical calculations show that for 8 > 2, i.e. depth of layer greater 





456 A. L, FLORENCE 


than the diameter of the disk, the contact stress 7,(r,0) is within | per 
cent of its value corresponding to a semi-infinite elastic medium (5c). In 
fact, for 8B = 2, 


7,4(7,0) = — * cpl 1-0099(1 —p2) ~§-+-0-0012(1 —p?)! +0-0001(1—p?)!]. 
T 


Again, for 8 = 2, 7’ = 5-387ya0a* compared with 7' = 5-333uaa* for 
B =o. The above expression for the contact stress when 8 = 2 was 
"730 


4pua 
-s/o, 


8-0 
6-0 
4+0F 
2:0) 














0 020406 08 1 12 14 16 18 20 22 
- 
a 


Fic. 4 


checked using the method of C. J. Tranter (8) in which a series representa- 
tion is employed instead of (6). A few curves have been plotted in Fig. 3 
to show the variation of the contact stress with layer depth. Those corre- 
sponding to B = } and } have been included to indicate the trend (the 
error for f(z) having been estimated at a value less than 10 per cent. 
Fig. 4 shows the variation of the torque with layer depth. 


3. Moment applied to the disk 
The vector for the moment shown in Fig. 5 lies in the x-axis directed 
away from the origin. 
Introduce Boussinesq—Papcovich functions ®, and ®, to describe the 
displacement field in the form (9) 
2uv = —grad ®,—zgrad®, +(3—4v)®, Z, (9) 
with V0, = V*O, = 0. 


The displacement vector v has components (u,, %g, u,) in cylindrical coordi- 
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nates and vy is Poisson’s ratio for the material and 2 is a unit vector in 
the direction of the z-axis. 

Let ®,(r,0,z) = d(r,z)sin@ and ®,(r,é,z) = ¢,(r,z)sin@ (10) 
Then from (9) and (10), using the strain-displacement and stress-strain 





a r 


elasti 


rigid 





relations, the following formulae for those components of displacement 
and stress concerned in the boundary conditions can be derived 


2 
2uu, 


me | 2 +(3—4»)b—4,sin 6 
fora—v) Pt 9 Th 
| ez éz® = az | 
a s|--% +(1—2»)$,—de sin # 


T= 5 —2% +(1—2v)dq, - d,|008 6 


a, = sin 0 


9 


Oz 


ve ~ 
~ 





in which by = Obp/ Ez. 
The assumed boundary conditions of the problem are 
u(r,6,0) = arsin@ (0 <r <a), o{r,0,0) = 0 (r> a), 

7, (r, 9,0) = r(r, 0,0) = 0, 

TAT, 0,h) = r.9(r,0,h) = ufr,0,h) = 0, (12) 
where « is the small rotation of the disk about its diameter. Unlike the 
torsion problem, it is assumed in (12) that no shear stresses exist on the 
contact surfaces of the elastic layer. Using formulae (11) the boundary 
conditions (12) for the functions ¢, and ¢, become 


(3—4v)4¢,—¢, =e ige 0,0 <r 2 @), 
£ (2(1—»)$,-43} = 0 (2 = 0,7 >a), 


(1—2v)d, — 2= 0 (z= 0), 


¢,=9 and ner +4, =O (z= 


~ 
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Since the functions ®, and ®, are harmonic, relations (10) and boundary 
conditions (13) lead to the following generalized form for functions ¢, 
and ¢,: 


wo 


a 2sinh fh : 
b= fA) ag Ep sinh Eh —2) (Er) dl, 
0 








= 2 sinh ¢h " 
$, = | 4@z55 —— 3gm cosh £(h—z)+- 
i) 


-+-(1—2v—£h coth €h)sinh £(h —z)}J, (Er) d€, 
where A(£) is determined by the dual integral equations 


oa 


| 4@ ogee ee = 0 <r <a), 
0 





€(2£h+-sinh 2£h) 1 


[ A@AEr)d—E = 0 (>a). (14) 
0 
Making the substitutions u = af, A(é) = {yaa/(1—v)}g(u), r = pa and 
h = Ba in equations (14) yields 


mw 


2sinh*Bu — 
—— mp (<p <1 
[ 9 Pre 28u+sinh 28u)° Jiloujen = p (0 <p <1) 


[ g(u)Jy(pu)du = 0 (p> 1). 


ri 


0 


Representation (6) as in the torsion problem leads to the Fredholm 
integral equation (7) for f(z), but now the kernel is 


K(p) =| v+sinhv.e~* cos po dv. 


v+sinh v cosh v 





0 
Adopting the procedure of § 2 for finding f(z) the contact stress o,(r, , 0) 
may be expressed in the form 


x 


o{r,0,0) = — oa [ (w)i(pu) dru sin @ 
cad 


M 
2"n! m! 
> ae — p2)n—t = > te 6. 
- v)? co (2n)!  (m— =myi amtt ex 


Also the moment mie to maintain a rotation « is 


M 
el ! m! 
a = eS 


~, (m—n)! 


+) Fam+1° 
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Some very elementary numerical calculations involving bounds show 
that when h = 3a the contact pressure o,(r,4,0) is between 1 per cent 
and 1} per cent greater than that corresponding to the semi-infinite 
medium. When B = o, 


3 
f(x) = z, o,(r, 0,0) = — - P sing, and M= Spoa™ 
3(1—v) 





m(1—v) (1—p?)! 
in agreement with G. N. Bycroft (2). 

In the boundary conditions (12) the assumption was made that the disk 
remains in contact with the elastic medium over the entire circle. This 
may be realized by a suitable superposition of the moment problem with 
that of a rigid disk pressed vertically into the elastic material. The latter 
case has been solved by N. N. Lebedev and Ia. 8S. Ufliand (10) and later 
by L. 1. Vorovich and Iu. A. Ustinov (11). 

The above analyses also apply to an infinite plate of thickness 2h with 
a rigid disk placed on either side of the plate exactly opposite each other 
and subjected to equal and opposite torques or moments. 

From the foregoing it can be concluded that for the analogous dynamic 
problems it is reasonable to assume that the contact stresses are those 
corresponding to the semi-infinite medium (with a periodic time factor) 
for stratum depths greater than twice the disk radius for torsional modes 
and three times the disk radius for bending modes. The corresponding 
stratum depths for the vertical mode are obtainable from the results of 
(10) and (11). 
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SUMMARY 

It is shown that, in general, Rayleigh waves in semi-infinite anisotropic media are 
either (i) dissipative, or (ii) not propagated at all. It is possible that non-dissipative 
Rayleigh waves can travel in certain discrete directions. These directions are 
probably the borders between the wedges in which the Rayleigh waves are dissipa- 
tive and the wedges in which they do not exist. 

If the free plane surface, however, is a plane of symmetry of the anisotropic 
material, Rayleigh waves, if they exist, are neither dissipative nor dispersive. The 
geometry of their propagation and the decay of their amplitudes is discussed in 
terms of the slowness and wave curves. 


1. Introduction and fundamental equations 
THE equations governing wave propagation in anisotropic elastic media 
were expressed in tensor form by Synge (1, 2), who also discussed the 
existence of Rayleigh surface waves in such media. He concluded that, 
in general, Rayleigh waves cannot be propagated, except, perhaps, in 
certain discrete directions. In this paper, a more general conclusion is 
obtained. It is shown that in anisotropic media Rayleigh waves, if they 
exist, are generally dissipative, except when the free plane boundary is 
also a plane of symmetry of the anisotropic material, when, if Rayleigh 
waves are possible, they are not dissipative. 

In a medium obeying Hooke’s law for sufficiently small displacements, 
the linear stress-strain relations in rectangular cartesian coordinates 2; are 


Tij = Cinq "pq? (1.1) 


where 7,; is the stress tensor, Eng the strain tensor, and the elastic coeffi- 


t 


cients ¢;,,, satisfy the symmetry relations 


c (1.2) 


tipa ~ “sipa = “tjqny = “pais? 
so that only twenty-one independent constants are involved. The suffixes 
can take the values 1, 2, or 3, and the summation convention for repeated 
suffixes is assumed. The strain tensor can be expressed in terms of the 
displacement vector u, by 

ty = (uw; 5+ U; 4). (1.3) 


Now at the Department of Applied Mathematics, The University of Sydney, N.S.W. 
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the commas indicating partial derivatives in the usual way. The equations 
of motion are 


otu : 
rs = Tyit+pX; (j = 1,2,3), (1.4) 


where p is the density, and X, is the body force per unit mass. Since only 
homogeneous media are considered, p and the ¢,,,, are constants, so that 
the three equations of motion obtained by combining these equations are 


a Cispq Mpg tPA. (1.5) 

We now use the plane wave approach of Synge (1) and Musgrave (3) 
to obtain the equation of the slowness surface. This method is simpler 
than the more powerful Fourier transform technique developed by Light- 
hill (4) and adopted by Buchwald (5, 6). Assume that there are no body 
forces, and that the displacement can be expressed in the form 


u; = A,exp[i(a,2,—wt)], ‘ (1.6) 
where A;, «,, and w are constants, w being the real circular frequency. 


If we substitute this expression in (1.5), we find that «,, must satisfy the 


equat‘on Q(a,) = det{M,,] = 0, (1.7) 


where Mj = Cijpq % %q—Pw*Ojy; (1.8) 


5,, being the Kronecker delta. The function Q(a,) is a homogeneous, sextic 
polynomial in a, w, and, from the relations (1.2), the matrix M;,, is 
symmetric. If the a’s are real, then the equation 2 = 0 represents 


the slowness surface for wave propagation in infinite media. 


2. The equation of Rayleigh waves 
Let the plane x, = 0 be the free boundary of a semi-infinite medium 
which occupies the space x, > 0. The boundary conditions are, therefore, 


T13 = Tes = Tgp = 0 (2.1) 
on #, = 0. Using (1.1), (1.2), and (1.3), these conditions can be expressed 
ns Cometng = 0 (= 1,3,3) (2.2) 
on %, = 0. 

We now look for solutions of (1.5) of the form (1.6) which satisfy the 
boundary conditions (2.2), and the amplitudes of which decay exponentially 
with increasing 25. If we assign given real values to «,, a, then a, can be 
obtained as the solution of the sextic equation Q(a,) = 0, each real root of 
which corresponds to an intersection of the straight line through the point 
(x,, %, 0), parallel to the a; axis, with the real surface Q = 0. Using the 
fact that the strain energy W = }c;;,, €;;€pq i8 always positive, it can be 
shown that Q = 0 is a closed surface. Hence there exists a curve T° in 
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the a,, a, plane, so that for values of «, and a, outside I’, the straight line 
parallel to the a, axis has no real intersections with 2 = 0, and the 
corresponding solution of the sextic Q(a,) = 0 consists of three pairs of 
complex conjugate roots for «,. Let af), where N takes the values 1, 2, 3, 
denote the three roots whose imaginary parts are positive. The three 
corresponding vectors Aj’) can now be found from the solution of the 


equation Ai) MW) = 0, (2.3) 


where the matrix M§*) is here determined by the given values of a, a, 
and the corresponding values of a‘). The rank of this matrix is generally 
2, so that the ratios Aj‘): AY”: AQ) are known in terms of its cofactors. 
The disturbance 3 

u; = - * 2a AN exp[i(al”)x, —wt)], (2.4) 


where f are arbitrary constants, satisfies the equation of motion, and 
decays exponentially as x, tends to infinity. Substituting this expression 
for u, in the boundary conditions (2.2), we obtain the three equations 


3 
> Cay fMASM a) = 0 (q = 1,2,3) 
N=1 


for the three unknowns f™. If 
K, 
then the equations (2.5) are consistent if 


D(o, «9, ) = det[K,x] = 0. 


om (N)(N) 
NW = Cg Aj”, 


This equation, homogeneous in «,, a, and w, was first obtained by Synge 
(1), and the above derivation is based on that given in his paper. 

For given real values of «,, a, the corresponding values of af) are 
complex, so that K,y is complex, and the equation D = 0 is an equation 
with complex coefficients, equivalent to two real equations 


D,(x, %,w) = 9; Dy{ox1, %g, v) = 9. (2.8) 
The two curves (2.8) in the «,, «, plane intersect only at a finite number 
of discrete points, and, therefore, except for some possible discrete direc- 


tions, the equation (2.7) has no solution when a,, a, are real. However, 
we may choose a, and a, to be complex and of the form 


a, = 7, +18), Mg == 2+ ty. (2.9) 


We now have two equations (2.7) for four unknowns. If we choose a 
direction in the z,, x, plane which makes an angle ¢ with the 2, axis, 
then for a plane wave travelling in this direction, 


f, = An, (2.10 a) 
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where A = tand. The relation 
8 = A, (2.10 b) 


must also hold, otherwise displacements on part of the wave front are 
exponentially large. This can be checked by making a rotation through 
an angle ¢ about the x, axis, when in the rotated coordinates the condition 
that amplitudes on a wave front z, = const. should not vary with 2, is 
given in (2.10b). 

If «,, a are rea! and outside [', the corresponding equation in a, has 
three roots whose imaginary parts are positive. We allow s, to increase 
from zero, and, up to a certain limit, three of the roots for a, have their 
imaginary parts positive. With s, above this limit the situation is not 
clear, as there may be only two roots with positive imaginary parts. 
However, there is a range in which s,, s, are positive, and the equation for 
x, has an equal number of roots above and below the real axis. We take s, and 
8, within this range and look for solutions of the equations (2.7) and (2.10). 
Such a solution may, or may not, exist for each given value of ¢. The above 
argument is only valid when ¢ is in the first quadrant, but it is always 
possible to arrange this by a suitable rotation of the z,, x, coordinates. 

Summarizing our conclusions: 

(i) there may be certain discrete values of ¢ for which the equation 
D(a, %,#) = 0 has a solution such that a, and a, are real; 

(ii) there may be ranges of values of ¢ in which this equation has solu- 
tions which are complex, corresponding to damped Rayleigh waves which 
decay exponentially with distance travelled; and 

(iii) for other ranges of ¢ the equation D = 0 may not have any solution 
at all. 

It seems likely that Case (i) is a borderline between Cases (ii) and (iii). 


3. Planes of symmetry 


THeorem. The equation D(a,,«.,w) = 01s purely real if the plane x, = 0 
is a plane of symmetry of the anisotropic material. 


It can be shown on the lines of (say) Sokolnikoff (7) that the plane 
x, = 0 is a plane of symmetry if and only if ¢,;,, = 0 when either one 
or three of the suffixes ¢, j, p, q take the value 3. There remain, at most, 
thirteen non-zero, independent constants, which are of the form C3333, C33,<, 
Carge, and C,,,,,, Where r, 8, t, wu can only take the values |, 2, and the sym- 


metry relations ¢4;,4 = Cjing = Cijgp = Cpqijs APE taken into account. From 


(1.8), May = Ciyrg jy — pes? (3.1) 


If Cj, is different from zero, then either i = g = 3, or i 4 3 and q ¥ 3, 
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so that M,, is a linear function of of. Similarly, M,,, M,,, M,,, and Mg, 
are also linear functions of of. On the other hand, 


Mis = Ciraq % %> (3.2) 
where now either i = 3 and g + 3, org = 3 andi + 3. Thus 


M,; = > ‘ (Cy133+Ca1ae)% %3- 


Similarly, the expressions for M,,, M,,, and M,, are also directly propor- 
tional to a3. It is easy to deduce that the equation 

Q(a;) = det{ M,,,] = 0, (3.3) 
is now a cubic in o§ with real coefficients. 

We choose a,, a, to be real and outside [, when the three values for 
a3 obtained from (3.3) are either (i) all real and negative, or (ii) one real 
and negative, and the other two complex and conjugate. From (2.3), 

AW) — A) A) 

we ME) My Me (ey Me)’ 
. a ; “S Fone : ~ 
ME MP) |My MEP) |My Mp! 

From (2.6), we have 
Kyy = Cysgi(ag AY) +0, Ay’), (3.5) 

with a similar expression for K,,, while 

Kyy = C3333 03. AZ") + C555, A}, (3.6) 
where NV = 1, 2, 3, summation over N is not implied, and the summation 
over ¢, t is for values 1, 2 only. We consider cases (i) and (ii) separately: 











(3.4) 


Case (i). [a§’}? is negative, and af’) purely imaginary. 

Using (3.1) and (3.2), it can be seen from (3.4) that A{’), Af’) can both 
be taken as real, with A$’) imaginary. Thus, from (3.5), both K,, and 
K,, are imaginary, and, from (3.6), Kj, is real. Hence the equation 

Ky, Ky Kys| 
Ky Keg Ks! 
is purely real. 

Case (ii). Without loss of generality, we may take | a" |* to be the real 
and negative root of (3.3), so that a} is purely imaginary. The other two 
roots [a{?’}* and [af]? are complex conjugates, whence a{” and af) have 
equal, positive, imaginary parts, while their real parts are equal in magni- 
tude, but opposite in sign. In other words, af? — —a{?’, the bars indicating, 
as usual, complex conjugates. 

As in case (i), Aj”, A} are real, and Af? is imaginary, so that K,,, Ky, 

5002 .56 Hh 
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are imaginary, and K,, is real. From (3.4), (3.1), and (3.2), we can deduce 
that A,, A, are even functions of a,, while A, is an odd function. Hence 
A® = A®, A® = A, a A® — a?) A®, and A? = —A®. Using these 
relations in (3.5) we deduce that K,, = —K,, and K,, = —K,., while, 
from (3.6), Ky, = Ay,. If we now multiply the first two rows of (3.7) by i, 
the first column consists only of real elements, while the elements of the 
second column are each the complex conjugate of the corresponding ele- 
ment in the third column. A brief examination shows that such a deter- 
minant takes only purely imaginary values, and, moreover, the elements 
in the second and third columns can be replaced by their real and imaginary 
parts, respectively. If we put L,, = iK,,, Ly, = 1Ky, Lg, = Kg, 
Lyg = re(t Ky), Legg = reli Kg), Lg. = re( Kg), Ly, = im(tK,,), 22, = im(iK,,), 
and L,, = im(K,,), then (3.7) reduces to the real equation 


det[ Ly] = 0. (3.7 a) 


This concludes the proof of the theorem. 

We have established that when z, = 0 is a plane of symmetry, the 
equation (2.7), D(a,,«,,w) = 0, has real coefficients and, therefore, repre- 
sents a curve in the real «,, a, plane. From (1.7), (1.8), (2.3), and (2.6) it 
follows that the function D is homogeneous in 4, a,, and w, so that we 
can obtain points on the curve by putting a, = Aa,, where A = tan¢, 
and then solve the resulting equation in a,/w. If there is such a solution, 
the corresponding vector (a,/w, «,/w) represeats in direction and inverse 
magnitude the phase velocity of Rayleigh waves, which are non-dissipa- 
tive. Since D is homogeneous in a, ag, and w, there is no dispersion. 
However, for some values of ¢, or ranges of values, the equation D = 0 
may not have a solution for real a, and a,, so that we would, in general, 
expect the curve given by this equation to consist of discontinuous 
segments. 

From (1.8) and (2.6) it follows that if a,, ag, a, satisfy (1.7), then —a,, 

%g, —a are also solutions of this equation. Since Q is here a function 
of o2, the roots corresponding to —a,, —a, are a{*’. Making a change of 
sign in a, and a, in all the equations in this section, it can be shown that 
if a,, a, satisiy D(a,,a.,w) = 0, then the equation D(—a«,, —a,,w) = 0 
is also satisfied. The curve D = 0 is, therefore, symmetric for reflections 
in the origin. 

The curve D = 0 can be regarded as the ‘slowness curve’ of the Rayleigh 
waves, in the sense that the slowness surfaces of Rayleigh waves are cylin- 
ders with generators parallel to the «, axis, which intersect the plane «, = 0 
on the curve D = 0. The wave curve (curve of equal phase at time ¢ = 1, 
consisting of waves from a source at the origin at ¢ = 0) is given in terms 
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of the parameters a,, a, by 


wD, wD, 


= (3.8) 


= —--—.-, Ze = asigeeenaerlceta » 

a D, +a D, ; a, D, +0 D, 
where a,, a, satisfy D(a,,a,,w) = 0, and D, = @D/@a,, D, = @D/éa,. As 
the waves are not dispersive, (3.8) can be written in the form 


a, = —D,/D,, z,= —D/D,, (3.9) 
where D,, = @D/éw. The wave curve given by (3.8) and (3.9) is the polar 
reciprocal of the slowness curve with respect to a certain circle. The form 
of the equation of the wave curve can be deduced either by a two-dimen- 
sional form of the method of Musgrave (3), or by the Fourier integral 
method due to Lighthill (4) as adapted by Buchwald (6). Using the latter 
technique, the following information is obtained. 

(i) To an observer at a point P having coordinates (2x,,7,,0), where 
r = (a?+-2)' is large, the waves from a source at the origin O are given 
by all N points S,, of the slowness curve at which the normal is parallel 
to OP. The corresponding vectors OS, represent the directions and in- 
verse velocities of the waves observed at P. 

(ii) The N intersections of the line OP with the wave curve represent 
the ‘group velocities’ of the waves observed at P, i.e. the velocities with 
which energy is propagated. Each of these intersections corresponds to 
one of the points S, of the slowness curve. 

(iii) The amplitudes of waves observed at P, corresponding to the point 
S,, of the slowness curve, vary inversely as r'. |grad D|. |k|!, where 


grad D| = (D}+-D})', 
and & is the curvature at S,,, given by 
ee Dy, Dy— 2D, D, D,+ Dy Dy (3.10) 
(Di +3)" 
provided that k + 0 at S,,. 

(iv) At a point of inflexion of the slowness curve, k = 0. It can be 
shown that the corresponding point of the wave curve is a cusp, and that 
in this case wave amplitudes vary inversely as |r|'. |\grad D|.|h|', where h 
is the coefficient of the third term of the Taylor expansion of the distance 
of a neighbouring point on the curve from the tangent at S,,. 

The proof and full discussion of the above results has been given in a 
previous paper (6), in which transversely isotropic media were investi- 
gated. In such a medium, the free plane is a plane of symmetry if and 
only if it is either parallel, or normal, to the axis, perpendicular to which 
there is symmetry of rotation. As particular examples of both these cases, 
the results of computations for a number of materials, as well as the graphs 
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of the more interesting slowness and wave curves, have been given in (6). 
Cusps were discovered in the wave curves of both zinc and apatite, and 
in some directions in these materials there are as many as three, or even 
five, Rayleigh waves. 

Stoneley (8) has investigated Rayleigh waves in particular directions 
of certain cuts of cubic crystals. However, such crystals contain several 
planes of symmetry, and for any of these planes, it is necessary to obtain 
both the slowness curve and the wave curve before the geometry of the 
Rayleigh waves can be properly described. 

Another problem in which the methods of this paper could be used is 
the case of a uniform magnetic field applied to a semi-infinite, isotropic, 
perfectly conducting, elastic solid. The anisotropy caused by the mag- 
netic field has been shown (9) to be rather similar to transverse isotropy, 
and there is the appropriate symmetry of reflection when the applied field 
is either normal or parallel to the free surface. Undamped Rayleigh waves 
can be propagated in either of these two cases, but if the magnetic field 
is inclined to the free surface, only dissipative Rayleigh waves are to be 
expected. 


Note added in proof 


Deresiewicz and Mindlin (10) have considered a particular case of waves 


on the surface of a monoclinic crystal. Using the notation of this paper, 
they took plane waves parallel to the x, axis in the case where the plane 
x, = 0, but not the free plane z, = 0, is a plane of symmetry of the 
material. Using the method of the section above, it can be shown that, 
with a, = 0, the equation (3.3) is here also a cubic in of. Moreover, it 
can be deduced that the corresponding Rayleigh wave equation 


D(O, x, w) = 0 


is purely real, both when the cubic equation possesses real, negative roots, 
and also when two of its roots are complex. However, the more general 
equation D(a,,«,,w) = 0 is here still compiex. 

We conclude that if the plane x, = 0 is a plane of symmetry, rather 
than the traction free plane xz, = 0, the direction parallel to the z, axis 
is a discrete direction in which undamped surface waves are possible, but 
that, in general, undamped surface waves in other directions cannot be 
propagated. It should be emphasized that roots of the equation 


D(0, xg, ) = 0 


must be sought in the region of complex values of 3, as well as real values, 
otherwise important solutions may be overlooked. 
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ON THE FLEXURE OF SOME PRISMATIC BEAMS 
By E. DEUTSCH 
(Professional School Klement Gottwald, Bucharest) 
[Received 12 December 1960} 


SUMMARY 
A method is given for solving the Saint-Venant flexure problem of a prismatic 
beam whose cross-section can be mapped conformally on a semicircle of unit radius. 
As an example the method is applied to the flexure problem of a prismatic beam 
whose cross-section is a semicircle of radius a. 


1. Introduction 
Amone the different methods for solving the Saint-Venant flexure problem 
of a prismatic beam, the method of Milne-Thomson (1) stands out for its 
remarkable elegance and applicability. In this, the shearing stress, dis- 
placement, torsion, and centre of flexure may be found directly in terms 
of a single flexure function, and the origin, as well as the direction of the 
coordinate axes in the plane of the cross-section, can be chosen arbitrarily. 

When the mapping function which maps conformally the cross-section 
of the beam on a unit circle is known, the method of Milne-Thomson also 
gives the expression for the flexure function. 

In this paper we shall show how the flexure function may be obtained 
for known mapping functions which map conformally the cross-section 
of the beam on a semicircle. 


2. Nomenclature 

O origin of coordinates taken in the plane of the fixed end. 

Ox, Oy, OZ three mutually perpendicular axes, the axis OZ being 
parallel to the generating lines of the beam. 
length of the beam. 
section of the free end Z = I. 
boundary of S. 
point in which OZ meets the plane Z = l, i.e. O'(0, 0, 1). 
centroid of S. 

tq = tate 


A {{ (y—yq)* dS, B= (f (x—z,,)* d8, 


8 


H J (x—2,)(y—yq) dS. 
W ; total load of a distribution of shear stresses applied over 


the plane Z = I. 


{Quart. Journ. Mech. and Applied Math., Vol. XIV, Pt. 4, 1961] 
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W,, W, 


y 


forces acting at O’ parallel to the z and y axes respec- 

tively—statically equivalent to the shear system W. 

M moment about OZ. 

P = W,+iW,. 

i ee 

a “soclemoe etd, A = AB—H’*. 

B = a,+ta, = a irs eer (bars are used to denote 
conjugate complex quantities). 
Young’s modulus. 

m shear modulus. 

v Poisson’s ratio. 

Tz, Fy, Os normal stress-com ponents. 

T xs) Tyg Tey Shearing stress-components. 

Les Ge coordinates of the centre of flexure. 


Ze = Lote 


Ob om 





3. Basic equations 
The non-vanishing stress-components are (1) 
Trg — tT yg = p[O(z) —p2—G?2*—122], (3.1) 
= E(Z—))[a,(x—2xg)+a,(y—Yq)], (3.2) 
where p = ir—}(1+v)(BZg+PezQ), (3.3) 
q=}(1+2r)8, r= dP, (3.4) 
7 being a real constant whose value is determined by (3.6) below and ®(z) 


is the so-called flexure function, holomorphic over the whole cross-section 
S and satisfying the boundary condition 


02 x,-, 02 b 0 RRO | yi 
(2) Ba) = (p-+-qi-+r2)8 = — (B+ Ge-+F2)2— (3.5) 


at the boundary C of S. 
The value of + is determined by the relation 
2M 
p» 


= }(K+RK)+2r1+3(1—2v)(BJ +B), 


I= JJ 2idS = —}i | 22* dz, 


Cc 


J = ee | basta sg ann 


K = 2 { | 2(z) dS = | 23@(z) dz. 
Ss 


of 
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We note that J is the polar moment of inertia with respect to O’. The 


centre of flexure is given by 

i 0 

z, == 2 J (3.10) 
where /, denotes the value of M which corresponds to 


t = }iv(B2g—Bz,). (3.11) 


Thus, once ®(z) is known, the flexure problem is completely solved. 


4. The flexure function 

We now show that condition (3.5) suffices to determine uniquely the 
holomorphic function ®(z) whenever the cross-section can be mapped 
conformally on a unit semicircle. Let 


z = w({) (4.1) 
be the function which maps conformally the region S on the unit semi- 
circle |{| < 1, y > Oin a {-plane ({ = £+%7). 

The function ®(z) can be expressed in terms of the variable {, so that 

D(z) = D[w(l)] = ¥(2), (4.2) 
where ¥(f) is holomorphic in the semicircle considered. 

Denote by y the semi-circumference |{| = 1, 7 >0 and by 8 the 
diameter » = 0, —1 < & < 1, y+ being the outline of the semicircle 
taken in positive sense. 

Let o = e be a point on y and € a point on 8. We have for the points 


caved z=w(s), = Ao), 


— Fao) do 
Cc 


ds’ 


Oz a) oz 
— = w (ce) — —= 
ds’ és 


and thus, from (3.5) we find 


(a )eo"(o) + 5 Ko-)a'(o) = &(0~)w'(o)[p+9qa(o-) +ra(o)}+ 


+ 0(0)'(o-[P-+Gu(o)+FH(o)}. (4.3) 


where accents denote differentiation with respect to the argument. 
For the points on 5 we have 


z= w(€), z= a), 


dé orn 
a = ws. a= #8) T, 


and thus from (3.5) we find 
b(E)eo' (E)—P(E)a'(E) = B(E)w'(E)[p+-qa€) +ra(E)]— 
—wE)ao'(E)[P+-qu(E)+Fa(E)]. (4.4) 
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hs: forces acting at O’ parallel to the z and y axes respec- 
tively—statically equivalent to the shear system W. 
M moment about OZ. 
P = W,+4W,. 
_ AW,—HW, . _ BW,—HW,. — AB—H? 
meee. aes, a= EA P A= AB #H?. 


B = a,+ta, = ts ee (bars are used to denote 





conjugate complex quantities). 
E Young’s modulus. 
shear modulus. 
Poisson’s ratio. 
on Gye M normal stress-components. 
zo Ty Tey Shearing stress-components. 
se. coordinates of the centre of flexure. 
eo Te+tY,- 
3. Basic equations 
The non-vanishing stress-components are (1) 
Tg — yg = p[O(z)—p2—gz*—r27], (3.1) 
o, => E(Z—I))[«,(x—2g)+a,(y—Ya)], (3.2) 
where p = ir—}(1+-v)(BiZg+BzQ), (3.3) 
q=}(1+2»)8, r= 38, (3.4) 
7 being a real constant whose value is determined by (3.6) below and ®(z) 


is the so-called flexure function, holomorphic over the whole cross-section 
S and satisfying the boundary condition 


ez -, 0% i. MRS SR pi 
(2) = — Ble) — = (pt+-git+ r2)8= — (P+ Ge + F2)2— (3.5) 
at the boundary C of 8S. 


The value of 7 is determined by the relation 


= = }(K+K)+2r1+3,(1—2r)(pJ +B), 


I= a) zi dS = —}i | 28* de, 


‘és 


J = ti JJ 2227 dS - dz, 


K = 2 [ { 2) dS = [ =20(2) az. 
Ss 


Cc 
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We note that J is the polar moment of inertia with respect to O’. The 
centre of flexure is given by 


6 
z, = &2—, l 
z, = 2 aP (3.10) 
where M, denotes the value of M which corresponds to 


+ = }ir(BE,—fe,). (3.11) 


Thus, once ®(z) is known, the flexure problem is completely solved. 


4. The flexure function 
We now show that condition (3.5) suffices to determine uniquely the 


holomorphic function ®(z) whenever the cross-section can be mapped 
conformally on a unit semicircle. Let 


z = w(¢) (4.1) 
be the function which maps conformally the region S on the unit semi- 
circle \{| < 1, n > O in a {-plane ({ = +77). 

The function ®(z) can be expressed in terms of the variable {, so that 
D(z) = D[w(L)] = ¥(2), (4.2) 
where y(f) is holomorphic in the semicircle considered. 
Denote by y the semi-circumference |{| = 1, 7 >0 and by 68 the 
diameter » = 0, —1 < & < 1, y+8 being the outline of the semicircle 
taken in positive sense. 


Let o = e® be a point on y and é a point on 5. We have for the points 


warned z=w(c), Z=a(o-), 


= i @ (a~*) do 


o ds’ 
and thus, from (3.5) we find 


Wo)e' (0) + Moa'(o?) = H(o-\w'(0)[p+ga(o-") +r(0)]-+ 


1 a ne -_ 
+= 0(0)a'(o)[P+Gu(o)+Fa(o-)]. (4.3) 
where accents denote differentiation with respect to the argument. 
For the points on 8 we have 
z=w(f), 2= a(€), 


eae 


and thus from (3.5) we find 
(Een (E)—P(E)a'(E) = G(E)w'(E)[p+-qalé)+ro(E)]— 
—wE)ao'(E)[P+-qe(E)+Fa(E)]. (4.4) 
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Putting ¥(l) = Wh(L)w'(2), 
the boundary conditions (4.3) and (4.4) may be written 
F(C)+'¥(C) = Ule) ony, 
¥()—¥(O) = Vig) ond, 
where 
U(f) = Ca(l ew’ (2)[p+-qa(l-*) + rw()]+- 
+£-ten(0)ca'(C-*)[P+-Gu(l) +Fa(L-)}, (4.8) 

V(C) = Led(C)ea’(L)[ p-+-gea(C) +-re(L)] —Leo(L) a’ (L)[P+-Geo(L) +-Fa(2)). 

(4.9) 
Clearly, we have a(t) = w(£). 

Oar flexure problem has thus been reduced to the task of ety a 
function ‘¥({), holomorphic in and on the semicircle |{| < 1, 7 > 0 whose 
real part assumes the values $U(c) on the semi- -circumference y and 
whose imaginary part assumes the values (1/2i)V(é) on the diameter 5; 
that is, we have been led to a Volterra problem (2) and, as is known (3), 
its solution is given by 


1 
1—@? 1 ] 4 
¥(0) = 5 [ v6, aetna | MTR) - 
6 -1 


o?+-1 
20” 





Since do = icd®, cos 6 = 
using (4.5), we find 


Hue) = 5 | Ue—;—— 5) do+ 
; 


aes a 


tf (O( ay ei) a, (4.10) 


-1 
which determines (¢); then from (4.2) we obtain ®(z). From (3.7), (3.8), 
and (3.9) we have 


I=-}i | w(o)a%(a-)e'(o) do—}i f mean f)w'(€)dé, (4.11) 
y ws 
~ 1 
J = | w(0)a%(o-)w'(o) do+ [ w(€)ea*(e)eu'(é) dé, (4.12) 


Y =~} 


1 
K = { w(o)a(o-)Yo)w'(o)do+ { w(€)B(EW(E\w"(E)4g. (4.13) 
-1 


Y 
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Thus relations (4.10), (4.11), (4.12), (4.13) constitute the solution of our 
flexure problem. 

If the — function w(Z) can be continued analytically into the 
semicircle |{| < 1, » < 0, formulae (4.10)-(4.13) can be put in a form 
which is ae more convenient for computations. In fact, from (4.9) 
we see that in this case V({) is holomorphic in the interior of y+. There- 
fore, making use of Cauchy’s formula and Cauchy’s theorem, we have 


“= fr "Ole 1) tin 7) Viol 3-sxga) # = VO. 
1 Y 


and thus (4.10) can be written 


CH(C)w'O) = VO+5- | (U(@)— Vie —>—=—73) do. (4.14) 
Y 


Similarly, taking into account the fact that the functions w({)@*({)w'(¢), 
w(L)d°(L)ew'(L), w(L)d(L)b(O)w"(L) are holomorphic in the interior of y+, 
the line nha in (4.11), (4.12), (4.13) along the diameter 5 can be 
replaced by line integrals along y. Thus we obtain 

I= -} ( w(a)w’(a)[a*(a0~!) —G*(a)] do, (4.15) 
Y 


J = | wl(c)w'(c)[a*(o-!) —d*(e)] do, (4.16) 


| 
y 
={ w(o)w' (a)ys(o)[G(o-!) — &(o)] do. (4.17) 
Y 


5. Torsion 
In the case of torsion we have 
W, = 0, W, = 0, M # 0, 
and thus 8 = 0. Now, from (3.3) and (3.4) it follows that 
p = fr, q=r= 0. (5.1) 
The function ¥({) is given by (4.10) or (4.14) where now 
U(E) = tr[ Ca(E-*)ea"(L) — Eh L)e0"(C ‘)], (5.2) 
V(C) = tr[Cei(C)w’(L) + Cen( L)ea"(C)). (5.3) 
The relation between ®(z) and the complex torsion function F(z) is 
(*, 5) (2) = rF (2), (5.4) 


which may be easily obtained by comparing the expressions for the 
shearing stresses, expressed in terms of ®(z) and F(z) respectively. 
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6. Example 


As an illustration of the application of the method described above, 
let us consider a beam whose cross-section is a semicircle of radius a. 
The mapping function is clearly 


z= at. (6.1) 
A simple computation gives 


(6.3) 
(6.4) 
(6.5) 


pair, g= 41428, r= HB. (6.6) 
pete 
From (6.1) we find 
w(e)= ale) = ac, ae-) = =, w'(e) = &'(o~) = a'(e) = a. 


(6.7) 
so that, from (4.8) and (4.9), we ret 


a(q+-7) 


U(c) = at|p+ p+ 3 “)4-a¢+ne|, 


V(c) = a*{2iro?+-a(q+r—gq—7)o*], 
since, from (3.3), p—p = 2ir. 
To obtain the flexure function, we use the relations 


J oe By Bet 
J aaa) oo 
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where in(;+4) - o(f+o+e +.. }} 
Now, from (4.14), we have 


¥(f) = ane +Setetito +t — 7) tart + 





+t : - 7) 4 eter) oy 


tn 2 


+ -a)--) e e—al 
(6.10) 


In terms of z the flexure function is 
2 * 2 - o P 
(2) = OPEB), Cara) 4 Pe 


— a 4 
+g q =(2—-§ 


in 


g4+r—q— a? zs 
; fe G+r—4 =A te -4-F/,2_@\|in aa 
str z 7 2 z a—z 
(6.11) 
To determine +r from (3.6) we first calculate the value of K from (4.17). 
Using (6.7) we find s 
K=¢ (1—o*)p(o) do, 
where y(c) is given by (6.10). 


Using the relations . 
jo%1da=0 (n #0), 


= in (n = 0), 
2 


oo" do = ———., 

. 2n+1 
Y 

where n is any positive or sae integer, and 


[a In n( 2) do st 2 a n(-*2) ion t, 
l—o — 


Y 


| “in(7 ad do = 2in, “f in - o 
Y Y 
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one finds, after some algebra, 


K= p+ Pin, ate(5— *) sa 
ris 





15 


2 


and from this we deduce 


K+K _ atr(s ah * 


2 
Substituting from (6.6) in (6.12) we get 


2a® te 
—-> (g+r+9+7). 


fe 3+ 2v)(B+8). 


From (3.6) we now find 
M atn(5 ai *) __ a u(3-+4v)(B+-B) 
208 15 
2W, _ 16W, 
EB a‘xK’ 
*)— 16apW,(3-+4y) 


or, since B+p = 


ok, Pat Bi (6.13) 


os li5rE 


M atru( ne 


which determines the value of 7. 
To find the centre of flexure, we insert in (6.13) the value of + given by 


3.11), ly, e 
{ ) namely, 2av ¢p : B) 32W, 
‘= —— = >... 
3a 3a°72 EF 
__ 16apW, 
157° E 


Putting W, = 3(P+ P) and using (3.10), one finds 


We find M, 2 (32+ »(40—n?)]. 


l6iapy 
— Feng (3n?+(40—a4)]; 
3n?-+1(40—n? 
whence x 0, Y,. ea 5 i (6.14) 
: L52r9(1 +-v) 
which agrees with the calculation of Stevenson (6), Sokolnikoff (4), Uflyand 
(7). 


In the case of torsion, putting g = r = 0, p = ir, we have from (6.11) 


rn 2 4) 
(2) — =a"(1 “) ie ’ ( 
= - 


which—taking into account relation (5.4)—agrees with the solution known 
for the complex torsion function (4, 8). 





ON THE FLEXURE OF SOME PRISMATIC BEAMS 
Finally, from (6.13) the torsional rigidity is found to be 


p=. at»(5—-) 


T 2 @ 


also given by Stevenson (6). 
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